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EXECUTIVE  SUMMARY 

Since  the  start  of  the  project  on  August  1,  1991,  all  of  the  tasks  described  in  the 
original  proposal  (see  Table  I)  have  been  addressed  successfully.  However,  due  to  the  fact 
that  years  2  and  3  were  funded  at  a  strongly  reduced  level.  Tasks  4,  5  and  6  could  only  be 
carried  out  at  a  feasibility  demonstration  level.  Without  the  ability  to  complete  the  full 
experimental  matrix  for  these  tasks,  several  issues  regarding  reproducibility,  specificity  and,  in 
particular,  applicability  (viz.  type  of  microorganisms  or  other  aerosol  components  that  can  be 
analyzed  with  this  technique)  remain  unanswered. 
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• 
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(Biocompound  Analysis) 
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X 
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Task  7; 

(Reporting) 

• 

X 

1 
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H 

X 

1 

1 
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In  order  to  document  the  experiments  carried  out  thus  as  completely  as  possible  and 
lay  the  groundwork  for  future  continuation,  instrumental  design  considerations  as  well  as 


underlying  theoretical  aspects  and  actual  experimental  observations  have  been  discussed  in  a 
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series  of  articles  and  conference  presentations  (see  Appendix  1).  Besides  confirming  the 
overall  feasibility  of  laser  pyrolysis/ionization  mass  spectrometry  of  single  organic  aerosol 
particles  by  means  of  ion  trap  methods,  the  results  of  the  work  indicate  the  potential  for 
several  spin-off  technologies,  including  precise  mass  determination  of  megadalton  ions  by 
means  of  nondestructive  optical  techniques  as  well  as  the  possibility  of  introducing  and 
trapping  large  ions  produced  by  atmospheric  pressure  ionization  methods,  e.g.,  electrospray 
ionization.  The  latter  topic  is  currently  the  focus  of  special  attention  as  it  appears  to  offer  an 
opportunity  for  continued  funding  under  ERDEC  auspices  while  utilizing  much  of  the 
expertise  and  instrumentation  developed  under  the  project  described  in  the  present  report. 
Finally,  it  should  be  mentioned  here  that  co-investigator,  Dr.  George  Hars,  who  returned  to 
his  native  country  (Hungary)  in  1993,  has  submitted  a  proposal  to  the  Joint  US/Hungarian 
fund  to  support  the  continuation  of  several  aspects  of  this  work  at  the  Department  of  Atomic 
Physics,  Technical  University  of  Budapest  in  direct  collaboration  with  the  University  of  Utah 
Center  for  Micro  Analysis  and  Reaction  Chemistry. 

From  the  work  reported  here,  the  following  conclusions  and  recommendations  can  be 
formulated. 

Conclusions: 

1)  Micron-  and  submicron-sized  particles  can  be  trapped  from  electrospray  sources  from 
vaporized  suspensions  in  H2O,  or  directly  fi'om  air. 

2)  Trapped  particles  can  be  pumped  down  to  high  vacuum  (<  1 0"^  torr)  conditions  without 
loss  of  stability. 

3)  Under  high  vacuum  conditions  particles  describe  a  star-like  pattern  in  the  equatorial 
plane  that  can  be  used  for  precise  mass/charge  ratio  determinations. 
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4)  Use  of  a  Nd  YAG  laser  beam  at  1 .06  ^lm  produces  intense  ion  signals  from  single 
styrene  particles  or  B.  subtilis  spores  which  appear  to  represent  complex  organic 
compounds. 

Recommendations: 

1)  Further  experiments  are  needed  to  optimize  Nd  YAG  laser  pulse  intensities,  beam 
alignment  and  choice  of  wavelengths  (the  1.06  pm  wavelength  can  be  doubled,  tripled 
or  quadrupled),  especially  in  order  to  increase  shot-to-shot  reproducibility. 

2)  Use  of  tandem  MS  techniques  should  enable  identification  of  the  presently  unknown 
ion  signals. 

3)  A  systematic  analysis  of  different  types  of  aerosol  particles,  supported  by  pattern 
recognition  methods  will  be  needed  to  evaluate  the  applicability  of  the  new  technique 
to  various  biodetection  problems. 


mXRODUCTION 


The  development  of  a  multifunction  (electrodynamic  balance  +  ion  trap  mass 
spectrometer)  Paul  trap  capable  of  trapping,  stabilizing  and  analyzing  microparticles  (ranging 
in  size  from  0. 1  to  several  tens  of  microns)  as  well  as  ions  (ranging  in  size  from  a  few  Daltons 
to  gigaDaltons)  was  undertaken  for  the  primary  purpose  of  obtaining  characteristic  mass 
spectra  of  single  microorganisms  and  other  respirable  aerosol  particles.  As  reported  at  the 
1992  ERDEC  Scientific  Conference^  successful  trapping  and  visualization  of  micron-,  or  even 
submicron-,  sized  particles  injected  by  means  of  various  aerosolization  methods  was 
accomplished  in  a  modified  Finnigan  MAT  type  ion  trap  electrode  assembly.  Optimization  of 
the  trapping  potential,  using  a  square  wave  drive  frequency,  was  shown  to  enable  evacuation 
of  the  chamber  to  less  than  10'^  torr  without  destabilizing  the  trapped  particle.  However, 
preliminary  attempts  to  produce  characteristic  mass  spectra  by  means  of  hitting  the  stabilized 
particle  with  the  focussed  beam  of  a  small  TEA  COj  laser  failed  to  produce  adequate  ion 
signal.  Serendipitously,  particles  trapped  at  pressures  <10"'*  torr  were  observed  to  exhibit  a 
previously  unreported  star-shaped  trajectory  in  the  equatorial  plane.  This  trajectory  was 
found  to  be  a  function  of  m/z  value. 

The  present  report  describes  the  evolution  of  a  novel,  nondestructive  optical  method 
for  determining  the  mass  of  a  single  megadalton  ion  (or  microparticle)  with  a  potentially  high 
degree  of  accuracy.  Secondly,  the  authors  report  the  installation  of  a  more  powerful  (Nd, 
YAG  type)  laser  and  the  results  of  preliminary  laser  fragmentation/ionization  (LFI)  MS 
studies  on  polystyrene  microparticles  and  Bacillus  subtilis  spores. 
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EXPERIMENTAL 


Instrumentation 

Figure  1  shows  the  experimental  setup.  A  modified  Finnigan  ITD  700  ion  trap 
electrode  assembly  is  used.  The  ring  electrode  has  three  radial  drill  holes  (3  mm  dia).  Two 
are  directly  in  line  with  each  other,  the  third  one  is  positioned  at  a  90  degree  angle.  The  two 
aligned  holes  carry  the  HeNe  laser  beam  through  the  center  of  the  trap.  The  perpendicular 
hole  is  used  to  detect  and  observe  particles  in  the  center  of  the  trap  levitated  by  means  of  a 
long  working  distance  optical  microscope  (Cohu  312985)  with  highly  sensitive  CCD  camera 
(Cohu  6410).  The  camera  can  also  be  operated  from  a  vertical  position  for  observing  the 
center  of  the  trap  through  the  central  hole  in  the  upper  end  cap.  The  ring  electrode  is  driven 
with  a  rectangular  voltage  signal  ranging  from  500  Hz  to  3  kHz  in  frequency  and  up  to  2  kV 
in  amplitude.  The  end  cap  electrodes  are  driven  by  a  power  supply  providing  both  DC  and/or 
AC  modulation  in  common  or  in  opposite  phase  mode.  A  UV  lamp  placed  underneath  the 
trap  projects  its  light  through  the  central  hole  of  the  lower  end  cap,  thus  illuminating  the 
center  of  the  trap.  The  pulsed  beam  (5  ns,  450  mJ)  of  the  Nd-YAG  laser  (Continuum, 

Surelite  J-10)  is  focussed  to  the  center  of  the  ion  trap  without  illuminating  any  part  of  the 
electrode  system.  The  trap  is  contained  within  a  UHV  vacuum  housing  with  suitable  view 
ports.  The  necessary  pumping  is  provided  by  a  170 1/s  TPU  170  Balzers  turbo  pump+  rotary 
forepump  combination.  During  evacuation  of  the  trap,  the  pressure  is  monitored  by  means  of 
a  Baratron  gauge. 
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Figure  1 .  The  experimental  setup  of  the  particle  analyzer  equipment. 

Experimental  Methods 

First  the  particles  are  suspended  in  water  and  then  introduced  into  the  ion  trap  by 
means  of  a  syringe  through  the  central  hole  of  the  upper  end  cap.  Because  of  the  high  electric 
field  in  the  trap  most  of  the  particles  will  be  spontaneously  charged.  The  initial  value  for  the 
mass/charge  ratio  is  about  10^  a.m.u./unit  charge.  Generally  more  than  one  particle  is  trapped. 
Particles  near  the  center  of  the  trap  are  being  visualized  by  means  of  the  HeNe  target  beam 
and  a  highly  sensitive  video  microscopy  system.  Applying  opposite  mode  DC  on  the  end  caps 
and  gradually  increasing  the  drive  frequency  of  the  ring  results  in  ejecting  all  particles  except 
one.  This  particle  will  be  centered  in  the  trap.  Subsequently  the  vacuum  system  is  gradually 
pumped  down  to  <10'^  torr.  In  order  to  avoid  arcing  in  the  trap  during  evacuation,  a 
rectangular  voltage  signal  is  used  to  drive  the  ring  electrode.  The  rectangular  signal  provides 
the  highest  possible  trapping  potential  at  limited  (200  V)  voltage  amplitude  at  high  drag 
conditions.  Below  10"^  torr  the  particle  reveals  highly  characteristic  and  stable  oscillation 
patterns.  In  the  equatorial  plane  of  the  ion  trap  the  particle  performs  an  elliptical  motion.  At 
specific  voltages  and  frequencies  the  rapidly  circulating  particle  can  be  observed  as  a  quasi 
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standing  oscillation  which  resembles  a  star  pattern  (see  Figures  2  and  3).  By  counting  the 
"branches"  of  the  star  and  measuring  the  corresponding  voltage  and  frequency,  the 
mass/charge  ratio  of  the  particle  can  be  calculated.  Subsequent  UV  radiation  to  induce 
stepwise  electron  loss  events  enables  the  absolute  mass  of  a  single  macro-ion  in  the  10^-10*" 
Dalton  range  to  be  determined. 


Figure  2.  Eight  branch  star  pattern  of  a  1 .0  Figure  3,  Nine  branch  star  pattern  of  a  1 .0 

pm  diameter  polystyrene  particle  to  the  pm  diameter  B.  stibtilis  spore  particle  in  the 

equatorial  plane  of  the  ion  trap.  equatorial  plane  of  the  ion  trap. 


Chemical  characterization  is  performed  in  the  same  Paul  trap,  while  operating  in  the  Ion  Trap 
Mass  Spectrometer  (ITMS)  mode.  Moreover,  since  the  mass  determination  procedure  is 
nondestructive,  the  same  microparticle  or  macro-ion  can  be  analyzed  while  still  stabilized  in 
the  trap.  First  He  buffer  gas  is  introduced  into  the  vacuum  system  to  10“*  torr.  Because  of  the 
increased  drag,  the  star  pattern  of  the  particle  shrinks  until  the  particle  appears  as  a  single 
bright  dot  in  the  center  of  the  trap.  Subsequently,  the  trapped  particle  is  fragmented  and 
ionized  by  means  of  a  single,  5  ns,  450  mJ  pulse  from  a  Nd  YAG  laser  operating  at  1 .06  pm. 
Simultaneous  with  the  laser  pulse  the  electronic  drive  of  the  ring  electrode  is  switched  from 
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»/2 

Figure  4.  LFI  mass  spectrum  of  a  1  pm  dia.  polystyrene  particle  (or  a  small  cluster  of 
particles)  showing  dominant  ion  peak  at  m/z  104,  believed  to  represent  the  styrene  monomer 
unit.  The  pattern  shown  is  not  veiy  reproducible  in  that  repeat  shots  may  show  a  similar  ion 
cluster  at  m/z  values  an5rwhere  between  m/z  104  and  110. 


m  w 


vJz 

Figure  5.  LFI  mass  spectrum  of  z.  Bacillus  subtilis  (or  a  small  cluster  of  several  spores) 
showing  a  more  complex  pattern  than  Figure  3 .  A  similar  type  of  pattern  was  found  in  about 
40-50%  of  the  B.  subtilis  spore  analyses.  Note  frequent  peak  intervals  of  15,  17,  18  and  26 
amu,  suggesting  dominant  organic  ion  signals. 


RESULTS  AND  DISCUSSION 


Mass  Determination  of  the  Particle 

A  theoretical  discussion  of  the  novel  mass  determination  principle  introduced  here  is 
based  on  the  pseudo-potential  model  of  the  ion  trap^.  In  this  model  the  charged  particle  is 
considered  to  be  in  a  spring-like  linear  force  field,  although  the  spring  constants  D,  and  are 
unequal.  For  mathematical  purposes  the  origin  of  a  Cartesian  coordinate  system  is  put  into  the 
center  of  the  trap.  The  x  and  y  axis  are  in  the  equatorial  plane  and  the  electrode  system  is 
rotationally  symmetric  around  the  vertical  z  axis.  Since  x  and  y  are  equivalent  they  will  be  often 
referred  as  r  (standing  for  radial  direction).  The  vector  of  average  force  affecting  the  charged 
particle  is  F. 

F  -  -DJ^xi*yf)  -  Dzk  (1) 


Here  i,  j,  k  are  the  unit  vectors  of  the  coordinate  system  in  x,  y  and  z  directions  respectively. 
The  values  of  D,.  and  are  derivable  from  the  model  and  are  as  follows: 


(2) 


where  and  are  the  spring  constants,  m  and  Q  are  the  mass  and  charge  of  the  particle,  V 
and  Q  are  the  voltage  amplitude  and  the  angular  frequency  of  the  drive  and  finally  r^  is  the 


radius  of  the  ring  electrode. 

Apparently,  the  "spring"  in  z  direction  is  four  times  stronger  than  in  r  direction.  (D^  = 
4Dj)  The  equation  of  motion  can  readily  be  solved.  The  x  and  y  components  are  identical.  The 
motion  in  each  direction  is  a  harmonic  oscillation.  The  oscillation  frequency  in  x,y  and  z 


directions  are  and  co^  respectively; 
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(3) 


CO 


m 


m 


These  are  the  so-called  secular  frequencies.  Obviously  the  secular  frequency  in  z  direction  is 
twice  higher.  Substituting  (2)  into  (3)  the  values  of  the  secular  frequencies  result  as; 


QV 

2  m  rla 


=v/2. 


QV 


m  rla 


(4) 


So  far  the  average  force  affecting  the  particle  was  considered  to  be  a  function  of 
position  only.  In  reality  this  average  force  is  the  result  of  several  "pushes  and  pulls"  so  the 
trajectory  of  the  particle  shows  a  fine  structure  which  is  synchronous  with  the  drive  frequency. 
If  the  drive  frequency/secular  frequency  ratio  is  an  integer  then  the  fine  structure  on  the 
trajectory  will  seem  to  be  standing.  In  the  case  of  the  elliptical  motion  in  the  equatorial  plane 
the  fine  structure  due  to  drive  frequency  will  result  in  a  standing  star  pattern  provided  that  the 
star  pattern  condition  is  met; 


«  =  Ci 


(5) 


HL  ^  H  V 

Q"  ^  rl 


(6) 


where  n  is  an  integer.  From  (5)  the  mass/charge  ratio  of  the  particle  can  readily  be  expressed; 
The  value  of  n  can  be  obtained  by  counting  the  branches  of  the  star.  By  means  of  (6)  the 
mass/charge  ratio  can  be  determined  by  adjusting  a  standing  (mostly  eight  branch)  star  pattern 
and  by  measuring  the  voltage  and  the  frequency  of  the  corresponding  drive. 
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Illuminating  the  negatively  charged  particle  with  very  faint  UV  light  the  electrons  will  be 
removed  from  the  particle  one  by  one^.  Some  second  succession  can  be  readily  adjusted  by 
using  faint  enough  UV  light.  If  the  particle  to  be  analyzed  is  in  a  standing  star  pattern  the  loss 
of  electrons  will  cause  the  pattern  to  spin  again.  Then  the  pattern  can  be  stabilized  by  reducing 
the  drive  frequency.  The  mass/charge  ratios  before  and  after  the  electron  loss  can  be  expressed 
by  means  of  (6): 

m  ^  n  V  m  ^  nV 

2  (Z-l)e  ■  2  U) 

where  Z  is  the  ionization  level,  e  is  the  elemental  charge  1 .6  x  10'*®  as,  Qi  and  Q2  are  neighbor 
drive  frequencies.  The  absolute  mass  can  be  calculated  from  (7): 


71  n  V  e 
m=- - 

2  rliol-al) 


(8) 


Since  the  particle  carrier  at  least  several  hundred  unit  charges  the  measurement  can  be  repeated 
a  number  of  times,  thus  making  the  method  more  dependable. 

Laser  Fragmentation/Ionization  MS 

The  LFI  mass  spectrum  of  a  1.0  pm  polystyrene  particle  illustrated  in  Figure  4  is 
characterized  by  the  occurrence  of  mass  peaks  at  every  nominal  m/z  value  within  the  selected 
mass  range.  Furthermore,  the  mass  intensities  are  quite  high  (typically  lO’-lO®  "ion  counts"  full 
scale).  Finally,  the  spectrum  appears  to  be  dominated  by  organic  compound  signals  in  view  of 
the  frequent  occurrence  of  14,  15  or  28  amu  intervals  rather  than  the  m/z  23  or  39  intervals 
characteristic  of  sodium  and  potassium  ion  contributions  or  the  12  amu  intervals  associated 
with  carbon  cluster  formation.  Unfortunately,  detailed  chemical  interpretation  of  the  mass 
spectrum  in  Figure  3  is  presently  impossible  since  the  signals  do  not  directly  correspond  to  the 
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expected  styrene  and  styrene-like  decomposition  products.  The  most  likely  explanation  is  that 
the  high  initial  concentrations  of  ions  in  the  ion  trap  are  causing  significant  shifts  in  the  mass 
scale  calibration.  Similarly,  the  B.  subtilis  spore  spectrum  in  Figure  5  also  provides  high  ion 
intensities  combined  with  a  current  lack  of  biochemical  interpretability.  Complex  secondary 
reactions,  e.g.,  due  to  excessive  laser  intensities  and/or  an  unknown  mass  scale  calibration  shifts 
due  to  space  charge  phenomena  in  the  trap  may  well  be  responsible  for  these  problems.  Before 
any  detailed  conclusions  can  be  drawn,  however,  it  will  be  necessary  to  confirm  the  organic 
nature  and  information  content  of  the  ion  signals  observed.  Since  the  ion  trap  system  can  be 
operated  in  a  tandem  MS  mode  it  should  be  possible  to  produce  daughter  ion  spectra  and  thus 
identify  some  of  the  major  peak  signals  in  Figures  4  and  5. 
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ABSTRACT 

A  novel  technique  for  laser  mass  spectrometry  of  individual  particles,  e.g. 
microorganisms,  is  being  developed.  Present  paper  gives  a  detailed  discussion  on  the 
theoretical  and  experimental  aspects  of  trapping  a  submicron  size  charged  particles  pressures 
from  atmospheric  down  to  <10'^  torr.  The  ability  to  trap  particles  under  UHV  conditions  has 
provided  a  new  opportunity  to  study  "ion  trajectories"  as  they  perform  the  solution  of  Mathieu 
equation. 

Individual  microparticles,  mainly  microorganisms,  have  been  aerosolized  and  charged 
by  a  quasi-electrospray  technique.  A  Paul  type  three-dimensional  quadrupole  "trap"  was 
constructed  to  combine  the  properties  of  an  EDB  (Electro-Dynamic  Balance)^  capable  of 
capturing  and  stabilizing  micro-sized  particles,  with  those  of  an  ITMS  (Ion  Trap  Mass 
Spectrometer),  capable  of  trapping  and  mass  selectively  detecting  ionic  species  up  to  several 
thousand  amu^^.  A  TEA  COj  laser  (300  mJ  per  200  [jsec  pulse)  with  focussing  optics 
designed  to  produce  a  50-80  p-m  beam  waist  through  the  center  of  the  trap  is  used.  A  typical 
analysis  cycle  starts  with  the  trap  in  EDB  mode,  thereby  enabling  capture  and  stabilization  of 
one  or  more  particles,  followed  by  evacuation  of  the  trap  to  high  vacuum  (<10‘^  torr). 

Subsequently,  pyrolysis/desorption/ionization  of  the  trapped  particles  will  be  performed 
with  a  single  COj  laser  pulse  while  switching  the  trap  over  to  ITMS  mode,  thereby  enabling 
mass  selective  analysis  of  the  ionic  species  produced,  including  quasimolecular  (e.g., 
cationized'*)  species.  The  present  paper  presents  interim  results  on  ionization  and  stabilization 
of  the  individual  organic  and  biological  particles. 
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INTRODUCTION 


Previous  work  in  our  laboratory  has  established  the  feasibility  of  COj  laser  pyrolysis 
of  individual,  50-100  pm  dia.  microparticles  suspended  (levitaied)  in  the  center  of  an 
electrodynamic  balance  (EDB)  chamber  at  ambient  pressure,  followed  by  GC/MS  analysis  of 
the  pyrolysis  vapor  products  with  an  ion  trap  mass  spectrometer  (ITMS)  system  at  high 
vacuum®.  A  2  m  long,  180  pm  i.d.  transfer  line  GC  column  connecting  the  EDB  and  ITMS 
chambers  functions  both  as  a  pressure  reducer  and,  via  balli.Mic  heating,  as  a  temperature 
programmed  GC  column.  Both  the  EDB  and  the  HMS  cell  are  Paul-type,  RF-only  traps 
operating  in  different  experimental  regimes.  In  order  to  investigate  the  feasibility  of  analyzing 
single  particles  in  the  1-50  pm  range  while  eliminating  transmission  losses  between. the  two 
cells,  a  fully  integrated  chamber,  performing  both  EDB  and  ITMS  functions  was  constructed 
(see  Figure  1).  Major  experimental  challenges  are:  (1)  to  quickly  stabilize  the  particles  under 
high  vacuum  conditions;  (2)  to  pyrolyze  these  particles  with  a  short  COj  laser  pulse;  and  (3) 
to  efficiently  trap  ionized  laser  pyrolysis  products. 


Figure  1.  Typical  electrical  conditions  for  EDB  and  ITMS  trap  operation. 

EXPERIMENTAL 

A  schematic  of  the  particle  levitation  apparatus,  as  presently  conceived,  is  displayed  in 
Figure  2.  Sample  introduction  is  from  the  top  of  the  trap  via  electrospray  from  a  MeOH 
solution,  KeNe  laser  illumination  projects  into  the  page  while  the  COj  laser  path  will  be  into 
the  trap  from  the  opposite  side.  Particles  are  detected  via  scattering  through  a  3  mm  hole  in 
the  ring  electrode  using  a  Lieca  MonoZoom  7  lens  and  Cohu  6400  series  CCD  camera.  The 
vacumn  system  is  pumped  by  a  110  1/s  Balzer  turbo  molecular  pump. 

The  electrospray  injector  was  constructed  from  a  syringe  and  a  high  voltage  power 
supply.  A  26,5  ga  needle,  cut  off  and  polished,  formed  the  spray  tip  and  particles  were 
introduced  into  the  syringe  as  a  suspension  in  methanol  (typically  2.5  x  10'^  weight  %),  The 
needle  tip  was  connected  via  a  20  M  Q  current  limiting  resistor  to  the  high  voltage  (-3  to  +3 
KV)  power  supply  while  the  electrically  insulated  plunger  was  operated  manually  via  light 
tapping  to  eject  small  pulses  of  liquid  for  the  electrospray.  The  spray  is  accelerated  into  the 
trap  via  the  high  voltage  gradient.  The  selection  of  trapping  fields  and  needle  potential  are 
conelated  to  yield  effective  and  efficient  trapping  of  the  charged  particles.  Utilizing  a  specific 
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Figure  2.  Particle  levitation  apparatus. 

needle  size,  geometry  and  solvent,  the  voltage  is  adjusted  to  give  a  desired  initial  droplet  size* 
which  is  slightly  larger  than  the  particle  size  (so  that  the  mass  does  not  decrease  much  during 
evaporation).  The  charge  level  is  slightly  below  the  Rayleigh  limit  using  a  2.2  KV  needle 
voltage.  This  allows  nearly  10*  charges  on  a  10[xm  particle. 

This  is  an  extremely  high  charge  to  mass  ratio  (-lO"*  charges  per  Da)  which  gives  these 
particles  a  high  immunity  to  destabilization  by  air  drafts  within  the  system.  This  allows  the 
system,  at  relatively  low  trapping  voltages  (<100  V  rms),  to  be  pumped  from  atmosphere  to 
10'^  torr  in  less  than  5  min. 

Using  the  2.2  KV  needle  voltage,  droplets  are  initially  trapped  using  a  ring  voltage  of 
“300  V  p-p  at  300  Hz.  While  the  end  cap  DC  voltage  difference  is  ~1  V.  As  the  solvent 
evaporates  particles  become  visible  near  the  center  of  the  trap.  The  electrospray  is  performed 
at  atmospheric  pressure  both  because  of  fear  of  glow  discharge  and  to  utilize  the  high  damping 
values  to  decellerate  the  charged  particles  entering  the  trap.  Although  the  electrospray  droplet 
is  accelerated  by  large  electric  fields  at  the  needle  tip,  the  atmospheric  damping  yields  a  large 
number  of  trapped  particles  from  a  single  pulse  of  liquid  (a  few  pL). 

Below  10"^  torr  the  frequency  and  voltage  of  the  trap  is  increased  up  to  about  1  KHz 
and  1000  V  p-p  in  order  to  avoid  Mathieu  instability  of  the  particle  since  the  effects  of  drag 
on  the  stability  map  are  not  seen  at  these  pressures.  Pressures  below  10'^  have  been  obtained 
without  destabilizing  the  particle. 

RESULTS  AND  DISCUSSION 

To  trap  a  charged  particle  the  following  conditions  should  be  met: 

1)  the  trapping  potential  well  should  be  deep  enough  in  both  r  and  z  directions  to  stabilize 
the  incoming  particle.  Typically  the  axial  dimension  is  critical  due  to  gravitational 
considerations.  The  particles  will  be  trapped  if 
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where  e  is  the  charge  on  the  particle,  is  the  RF  amplitude,  is  the  axial  dimension,  Q  is 
the  angular  frequency  and  A  is  the  damping  time  constant.  Additionally  must  be 
sufficiently  large  to  accomodate  additional  particle  perturbations  due  to  electrical  noise 
molecular  collisions  at  low  pressures  and  due  to  air  currents  at  higher  pressures. 
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On  the  other  hand,  the  parameter  of  the  particle  may  not  exceed  the  Mathieu 
stability  threshold. 
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This  threshold  for  A=0  is  0^=0.909  and  is  increasing  with  the  second  power  of  the  drag 
parameter  D=A/Q.  Figure  3  illustrates  the  stability  map  for  two  values  of  the  drag  parameter. 


Figure  3.  Mathieu  equation  stability  map  for  two  values  of  the  damping  factor  D. 

3)  Sufficient  damping  of  injected  particles  must  be  present  in  order  to  insure  sufficient 
numbers  of  particles  are  trapped  in  the  potential  well.  Figure  4  indicates  the  damping 
time  constant  as  a  function  of  particle  size  and  pressure.  When  A  is  greater  than  the 
trapping  firequency  Q,  energy  loss  to  drag  is  sufficiently  rapid  that  high  capture 
efficiencies  may  be  expected. 
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Figure  4.  Particle  damping  factor 
as  a  function  of  gas  pressure. 
Computed  for  air  at  20  C. 


These  three  conditions  yield  competing  criteria  for  particle  stability.  Large  damping 
time  constants  are  required  to  decellerate  the  electrosprayed  particles,  but  they  lead  to  low- 
trapping  potentials  so  particles  are  easily  destabilized.  Large  voltages  and  low  frequencies 
yield  deep  potential  wells,  but  may  lead  to  particle  instability.  Ultimately  one  selects  a 
compromise  to  best  accomodate  a  particular  m/z  value  being  produced  by  the  electrospray. 
Clearly  the  highest  e/m  ratio  give  the  greatest  degree  of  particle  control.  During  sample 
introduction,  several  particles  may  be  trapped.  The  trapping  voltages  and  frequencies  are 
manipulated  to  obtain  a  single  particle  in  the  trap.  It  should  be  noted  that  while  multiple 
particles  may  be  simultaneously  stabilized,  electrostatic  repulsion  between  them  prevents  either 
particle  being  damped  to  the  null  point  of  the  trap. 

After  stabilizing  a  single  particle,  pump  dewn  is  made  gradually.  In  the  hundred  torr 
range,  speeds  up  to  about  2  torr/sec  pressure  reduction  will  not  to  blow  out  the  particle.  At 
about  20  tort  the  trapping  voltage  should  be  lowered  to  150  Vrms  to  avoid  arcing.  When  the 
pressure  goes  to  10'^  torr  range,  the  frequency  of  the  drive  is  increased  to  avoid  Mathieu 
instability,  while  the  trapping  voltage  is  also  increased  to  about  600  Vrms  to  maintain  sufficient 

trapping  potential.  ^ 

Figure  5  illustrates  a  single  bacillus  subtillus  cell  (~3  jxm  is  diameter)  suspended  at  10' 
torr.  The  particle  trajectory  forms  a  Lysijou  figure  due  to  the  superposition  of  different  axial 
and  radial  oscillation  frequencies.  The  amplitude  is  approximately  100  |xm  across. 

Figure  6  illustrates  the  same  particle  at  10'^  immediately  (<  5s)  after  pressuring  the 
chamber  to  lO'^  torr  and  pumping  to  back  out.  The  trajectory  of  Figure  4  will  be  fully 
developed  after  a  number  of  minutes.  In  this  case,  the  particle  has  been  mechanically  excited 
by  vibrating  the  electrode  assembly.  The  particle  location  is  localized  within  50  p.m,  which 
is  adequate  for  precise  laser  targeting. 
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Figure  5.  Bacillus  subtillus  cell  suspended 
at  10-7  torr.  Image  particle  motion  over  1 
video  frame  (1/30  s). 


Figure  6.  Same  particle  as  Figure  5, 
immediately  after  application  of  damping 
gas. 


REMAINING  CHALLENGES 


Laser  Pyrolysis/Desorption/Ionization 

A  pulsed  TEA  COj  laser  will  be  used  with  dual  pulse  time  option  (10  ms  and  200  ms) 
and  capable  of  producing  300  mJ  pulses  in  TEMqq  mode.  Previous  work  with  in-source  COj 
laser  pyrolysis  of  biological  samples  has  conclusively  established  the  occurrence  of 
quasimolecular  ion  formation,  particularly  cationized  species,  of  relatively  large  biomolecules 
due  to  laser  desorption/ionization  processes^ 

Mass  Spectrometric  Detection  of  Product  Ions 


This  task  utilizes  additional  high  frequency,  ion  trap  type  mass  scanning  electronics  and 
wil  require  that  EDB  electronics  be  shut  off  in  <10  ms  prior  to  the  laser  pulse  to  prevent 
destabilization  of  resulting  ions.  Considerable  optimization  to  obtain  sufficient  ion  yield 
without  excessive  space  charge  will  be  required. 

Detection  and  Identificaiion  of  Biological  Aerosol  Particles 

The  complexity  of  biological  materials  require  subsequent  analysis  which  will  focus  on 
several  areas:  (1)  the  usage  of  pattern  recognition  and  multivariate  data  analysis  methods  to 
extract  "finger  prints"  of  various  biological  materials.  (2)  The  usage  of  MS/MS  techniques  to 
distinguish  the  origin  of  various  ions  which  may  be  associated  with  particular  materials.  (3) 
The  expansion  of  the  mass  range  to  higher  m/z  values  to  enable  examination  of  particular 

molecular  ions  which  in  themselves  may  be  more  highly  characteristic  of  various  biological 
samples. 
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CONCLUSIONS 


It  is  believed  that  10'^  torr  will  be  a  sufficiently  low  system  pressure  to  perform  the 
laser  pytolysis  and  mass  spectrometry  on  a  particle.  If  lower  pressures  are  required, 
stabilization  of  various  particles,  including  1 3  pm  paper  mulberry  pollen,  latex  microspheres 
below  .93  M-m  and  bacillus  subtillus  organisms  and  spores,  have  been  obtained  at  pressures 
down  to  10"^  torr,  but  achieving  these  pressures  may  require  special  measures  to  contain  the 
size  of  the  particle  trajectory  and  to  remove  condensate  (primarily  HjO)  in  the  vacuum  system 
that  limits  pumpdown  speed  below  10’^  torr. 

The  expansion  of  particle  trajectories  below  10’^  torr  is  due  primarily  to  thermal  and 
electrical  noise  which  is  no  longer  damped  out  by  the  gas  damping  forces.  In  this  region  the 
"Matthieu”  trajectories  are  readily  observed  (see  Figure  5).  It  is  not  clear  if  these  trajectories 
may  be  damped  by  electrical  means,  but  initial  efforts  indicated  that  if  damped  by  pressures 
on  the  order  of  10'^  torr,  the  particle  can  be  stabilized  back  at  the  center  of  the  trap.  Expansion 
of  the  trajectory  after  the  pressure  is  back  to  10’’  may  take  several  minutes,  so  that  gas  leaks 
provide  an  effective  stabilization  compatible  with  high  vacuum  pyrolysis  experiments. 

It  was  anticipated  that  the  position  would  be  less  localized  at  reduced  pressure  and  that 
additional  means  would  be  necessary  to  maintain  the  particle  at  the  trap  center  for  laser 
targeting.  It  was  not  anticipated  that  it  would  be  so  stable  to  such  a  low  pressure.  The  ability 
to  trap  particles  under  these  conditions  has  ‘provided  a  new  opportunity  to  study  "ion 
trajectories"  as  a  function  of  pressure  and  trapping  conditions.  Pursuit  of  this  area  could 
provide  valuable  insights  into  resolution  enhancements  provided  by  resonance  ejection’  and  into 
the  process  of  collision  induced  dissociation  in  Ion  Trap  MS  instruments. 

Finally,  it  was  not  anticipated  that  the  visibility  of  submicron  particles  will  be  so  good. 
This  is  mainly  due  to  the  high  sensitivity  of  the  CCD  camera  for  the  wavelength  of  HeNe  laser 
emission  (632  nm). 
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Abstract: 


From  the  Abstract  Book  for  the  American  Chemical  Society  Meeting, 
Chicago,  IL,  August  1993. 


CHEMICAL  CHARACTERIZATION  OF  SINGLE  AEROSOL  PARTICLES  VIA  COMBINED 
ELECTRODYNAMIC  BALANCE/ION  TRAP  MASS  SPECTROMETRY  TECHNIQUES.  Neil 
S.  Arnold,  George  Hars  and  Henk  L.C.  Meuzelaar;  Center  for  Micro  Analysis  and  Reaction  Chemistry, 
University  of  Utah,  Salt  Lake  City,  UT  84112. 

The  electrodynamic  balance  (EDB)  has  been  used  for  various  studies  of  individual  aerosol 
particles  with  respect  to  a  variety  of  physical  and  optical  properties.  Ion  Trap  mass  spectrometry 
requires  the  same  electrode  configuration  although  different  operating  voltages  are  used.  The  use  of 
combined  EDB  and  Ion  Trap  techniques  for  chemical  evaluation  of  single  aerosol  particles  allows  initial 
particle  evaluation  with  respect  to  size  and  weight  to  be  obtained  via  optical  techniques  during  initial 
particle  suspension,  while  laser  pyrolysis/desorption  is  used  to  generate  ionic  species  which  are  mass 
analyzed  for  chemical  characterization  of  the  particle  via  Ion  Trap  mass  spectrometry.  Optical  detection 
and  stabilization  of  individual  aerosol  particles  may  be  performed  for  particles  less  than  1  micron  to  over 
100  microns  in  diameter.  Novel  methods  for  determining  the  mass  of  these  particles  prior  to  laser 
pyrolysis/desorption  will  be  described  as  well  preliminary  mass  spectrometric  evaluation  of  the  products 
used  for  particle  characterization. 
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ABSTRACT 

A  multifunction  Paul  trap  is  described  capable  of  trapping,  stabilizing  and  analyzing 
electrostatically  charged  microparticles  as  well  as  a  broad  range  of  ions  and  macro-ions. 
Typically,  particles  in  the  0.1-10  p.m  range  are  introduced  by  aerosolization  from  an  aqueous 
suspension,  although  quasi-electrospray  and  dry  powder  type  introduction  methods  can  also  be 
used.  A  new  particle  trajectory  pattern,  observed  in  the  equatorial  plane,  was  found  to  offer  a 
nondestructive,  optical  method  for  determining  the  m/z  value  of  microparticles  and  macro-ions, 
with  a  present  accuracy  of  1:10^  and  a  potential  maximum  resolution  of  1:10®.  Recent  addition 
of  a  more  powerful  laser  (Nd  YAG,  operating  at  1.06  pm)  has  enabled  us  to  generate  laser 
fragmentation/ionization  mass  spectra  of  1  pm  dia  polystyrene  particles  and  of  Bacillus  subtilis 
spores.  Current  shot-to-shot  reproducibility  is  still  unsatisfactory  and  some  as  yet  unresolved, 
mass  calibration  problems  have  been  encountered.  Nonetheless,  the  high  intensity  and  apparent 
complex  organic  nature  of  the  ion  signals  obtained  might  herald  the  emergence  of  a  novel  MS 
technique  for  chemical  and  physical  characterization  of  single  microorganisms  and  other 
components  of  respirable  aerosols. 


INTRODUCTION 

The  development  of  a  multifunction  (electrodynamic  balance  +  ion  trap  mass 
spectrometer)  Paul  trap  capable  of  trapping,  stabilizing  and  analyzing  microparticles  (ranging 
in  size  from  0.1  to  several  tens  of  microns)  as  well  as  ions  (ranging  in  size  from  a  few  Daltons 
to  gigaDaltons)  was  undertaken  for  the  primary  purpose  of  obtaining  characteristic  mass  spectra 
of  single  microorganisms  and  other  respirable  aerosol  particles.  As  reported  at  the  1992 
ERDEC  Scientific  Conference^  successful  trapping  and  visualization  of  micron-,  or  even 


submicron-,  sized  particles  injected  by  means  of  various  aerosolization  methods  was 
accomplished  in  a  modified  Finnigan  MAT  type  ion  trap  electrode  assembly.  Optimization  of 
the  trapping  potential,  using  a  square  wave  drive  frequency,  was  shown  to  enable  evacuation 
of  the  chamber  to  less  than  10"^  ton  without  destabilizing  the  trapped  particle.  However, 
preliminary  attempts  to  produce  characteristic  mass  spectra  by  means  of  hitting  the  stabilized 
particle  with  the  focussed  beam  of  a  small  TEA  COj  laser  failed  to  produce  adequate  ion  signal. 
Serendipitously,  particles  trapped  at  pressures  <10'^  ton  were  observed  to  exhibit  a  previously 
unreported  star-shaped  trajectory  in  the  equatorial  plane.  This  trajectory  was  found  to  be  a 
function  of  m/z  value. 

The  present  report  describes  the  evolution  of  a  novel,  nondestructive  optical  method  for 
determining  the  mass  of  a  single  megadalton  ion  (or  microparticle)  with  a  potentially  high 
degree  of  accuracy.  Secondly,  the  authors  report  the  installation  of  a  more  powerful  (Nd,  YAG 
type)  laser  and  the  results  of  preliminary  laser  fragmentation/ionization  (LFI)  MS  studies  on 
polystyrene  microparticles  and  Bacillus  subtilis  spores. 

EXPERIMENTAL 


Instrumentation 

Figure  1  shows 
the  experimental  setup. 
A  modified  Finnigan 
ITD  700  ion  trap 
electrode  assembly  is 
used.  The  ring 
electrode  has  three 
radial  drill  holes  (3 
mm  dia).  Two  are 
directly  in  line  with 
each  other,  the  third 


Foreground  -  HeNe  lasw  ^  , 

Backoraund  -  Nd-YAG  laser  alectrospray 


Turtxj  Moiecuiar  Pumping  Roughing  +  Pressure  Control 


one  is  positioned  at  a  pjgure  1.  The  experimental  setup  of  the  particle  analyzer 
90  degree  angle.  The  equipment, 
two  aligned  holes 

carry  the  HeNe  laser  beam  through  the  center  of  the  trap.  The  perpendicular  hole  is  used  to 
detect  and  observe  particles  in  the  center  of  the  trap  levitated  by  means  of  a  long  working 
distance  optical  microscope  (Cohu  312985)  with  highly  sensitive  CCD  camera  (Cohu  6410). 
The  camera  can  also  be  operated  firom  a  vertical  position  for  observing  the  center  of  the  trap 
through  the  central  hole  in  the  upper  end  cap.  The  ring  electrode  is  driven  with  a  rectangular 
voltage  signal  ranging  from  500  Hz  to  3  kHz  in  frequency  and  up  to  2  kV  in  amplitude.  The 
end  cap  electrodes  are  driven  by  a  power  supply  providing  both  DC  and/or  AC  modulation  in 
common  or  in  opposite  phase  mode.  A  UV  lamp  placed  underneath  the  trap  projects  its  light 
through  the  central  hole  of  the  lower  end  cap,  thus  illuminating  the  center  of  the  trap.  The 
pulsed  beam  (5  ns,  450  mJ)  of  the  Nd-YAG  laser  (Continuum,  Surelite  J-10)  is  focussed  to  the 
center  of  the  ion  trap  without  illuminating  any  part  of  the  electrode  system.  The  trap  is 
contained  within  a  UHV  vacuum  housing  with  suitable  view  ports.  The  necessary  pumping  is 
provided  by  a  170  1/s  TPU  170  Balzers  turbo  pump+  rotary  forepump  combination.  During 
evacuation  of  the  trap,  the  pressure  is  monitored  by  means  of  a  Baratron  gauge. 


Experimental  Mt^thods 


First  the  particles  are  suspended  in  water  and  then  introduced  into  the  ion  trap  by  means 
of  a  syringe  through  the  central  hole  of  the  upper  end  cap.  Because  of  the  high  electric  field 
in  the  trap  most  of  the  particles  will  be  spontaneously  charged.  The  initial  value  for  the 
mass/charge  ratio  is  about  10“  a.m.u./unit  charge.  Generally  more  than  one  particle  is  trapped. 
Particles  near  the  center  of  the  trap  are  being  visualized  by  means  of  the  HeNe  target  beam  and 
a  highly  sensitive  video  microscopy  system.  Applying  opposite  mode  DC  on  the  end  caps  and 
gradually  increasing  the  drive  frequency  of  the  ring  results  in  ejecting  all  particles  except  one. 
This  particle  will  be  centered  in  the  trap.  Subsequently  the  vacuum  system  is  gradually  pumped 
down  to  <10"*  torr.  In  order  to  avoid  arcing  in  the  trap  during  evacuation,  a  rectangular  voltage 
signal  is  used  to  drive  the  ring  electrode.  The  rectangular  signal  provides  the  highest  possible 
trapping  potential  at  limited  (200  V)  voltage  amplitude  at  high  drag  conditions.  Below  lO"^  torr 
the  particle  reveals  highly  characteristic  and  stable  oscillation  patterns.  In  the  equatorial  plane 
of  the  ion  trap  the  particle  performs  an  elliptical  motion.  At  specific  voltages  and  frequencies 
the  rapidly  circulating  particle  can  be  observed  as  a  quasi  standing  oscillation  which  resembles 
a  star  pattern  (see  Figures  2  and  3).  By  counting  the  "branches"  of  the  star  and  measuring  the 
corresponding  voltage  and  frequency,  the  mass/charge  ratio  of  the  particle  can  be  calculated. 
Subsequent  UV  radiation  to  induce  stepwise  electron  loss  events  enables  the  absolute  mass  of 
a  single  macro-ion  in  the  10^-10^“  Dalton  range  to  be  determined. 


Figure  2.  Eight  branch  star  pattern  of  a  1.0  Figure  3.  Nine  branch  star  pattern  of  a  1.0 

pm  diameter  polystyrene  particle  to  the  pm  diameter  B.  subtilis  spore  particle  in  the 

equatorial  plane  of  the  ion  trap.  equatorial  plane  of  the  ion  trap. 


Chemical  characterization  is  performed  in  the  same  Paul  trap,  while  operating  in  the  Ion  Trap 
Mass  Spectrometer  (ITMS)  mode.  Moreover,  since  the  mass  determination  procedure  is 
nondestructive,  the  same  microparticle  or  macro-ion  can  be  analyzed  while  still  stabilized  in  the 


trap.  First  He  buffer  gas  is  introduced  into  the  vacuum  system  to  torr.  Because  of  the 
increased  drag,  the  star  pattern  of  the  particle  shrinks  until  the  particle  appears  as  a  single  bright 
dot  in  the  center  of  the  trap.  Subsequently,  the  trapped  particle  is  fragmented  and  ionized  by 
means  of  a  single,  5  ns,  450  mJ  pulse  from  a  Nd  YAG  laser  operating  at  1.06  jxm. 
Simultaneous  with  the  laser  pulse  the  electronic  drive  of  the  ring  electrode  is  switched  from  the 
mass  determination  to  the  ITMS  mode  in  order  to  record  a  mass  spectrum  of  the  ion  fragments 
created  by  the  laser  shot  (Figures  4  and  5). 

Figure  4.  LFI  mass  spectrum 
of  a  1  [xm  dia.  polystyrene 
particle  (or  a  small  cluster  of 
particles)  showing  dominant 
ion  peak  at  m/z  104,  believed 
to  represent  the  styrene 
monomer  unit.  The  pattern 
shown  is  not  very 
reproducible  in  that  repeat 
shots  may  show  a  similar  ion 
cluster  at  m/z  values 
anywhere  between  m/z  104 
and  110. 
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Figure  5.  LFI  mass  spectrum  of 
a  Bacillus  subtilis  (or  a  small 
cluster  of  several  spores) 
showing  a  more  complex  pattern 
than  Figure  3.  A  similar  type  of 
pattern  was  found  in  about  40- 
50%  of  the  B.  subtilis  spore 
analyses.  Note  frequent  peak 
intervals  of  15,  17,  18  and  26 
amu,  suggesting  dominant 
organic  ion  signals. 

RESULTS  AND  DISCUSSION 
Mass  Determination  of  the  Particle 

A  theoretical  discussion  of  the  novel  mass  determination  principle  introduced  here  is 
based  on  the  pseudo-potential  model  of  the  ion  trap^  In  this  model  the  charged  particle  is 
considered  to  be  in  a  spring-like  linear  force  field,  although  the  spring  constants  D,  and  are 
unequal.  For  mathematical  purposes  the  origin  of  a  Cartesian  coordinate  system  is  put  into  the 
center  of  the  trap.  The  x  and  y  axis  are  in  the  equatorial  plane  and  the  electrode  system  is 
rotationally  symmetric  around  the  vertical  z  axis.  Since  x  and  y  are  equivalent  they  will  be 
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often  referred  as  r  (standing  for  radial  direction).  The  vector  of  average  force  affecting  the 
charged  particle  is  F. 

F  =  -D/xi+yj)  -  (1) 


Here  i,  j,  k  are  the  unit  vectors  of  the  coordinate  system  in  x,  y  and  z  directions  respectively. 
The  values  of  and  D^.  are  derivable  from  the  model  and  are  as  follows: 
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where  Dj.  and  are  the  spring  constants,  m  and  Q  are  the  mass  and  charge  of  the  particle,  V 
and  Q  are  the  voltage  amplitude  and  the  angular  frequency  of  the  drive  and  finally  r,,  is  the 
radius  of  the  ring  electrode. 


Apparently,  the  "spring"  in  z  direction  is  four  times  stronger  than  in  r  direction.  (D^  = 
4Dj)  The  equation  of  motion  can  readily  be  solved.  The  x  and  y  components  are  identical.  The 
motion  in  each  direction  is  a  harmonic  oscillation.  The  oscillation  frequency  in  x,y  and  z 
directions  are  co^=cOy=cOj  and  co^  respectively: 
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These  are  the  so-called  secular  frequencies.  Obviously  the  secular  frequency  in  z  direction  is 
twice  higher.  Substituting  (2)  into  (3)  the  values  of  the  secular  frequencies  result  as: 
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So  far  the  average  force  affecting  the  particle  was  considered  to  be  a  function  of  position 
only.  In  reality  this  average  force  is  the  result  of  several  "pushes  and  pulls"  so  the  trajectory  .  _ 
of  the  particle  shows  a  fine  structure  which  is  synchronous  with  the  drive  frequency.  If  the  - 
drive  frequency/secular  frequency  ratio  is  an  integer  then  the  fine  structure  on  the  trajectory  will  -  - 
seem  to  be  standing.  In  the  case  of  the  elliptical  motion  in  the  equatorial  plane  the  fine  “ 
structure  due  to  drive  frequency  will  result  in  a  standing  star  pattern  provided  that  the  star 
pattern  condition  is  met:  _3 

I*  =  Q  (5) 

where  n  is  an  integer.  From  (5)  the  mass/charge  ratio  of  the  particle  can  readily  be  expressed: 

The  value  of  n  can  be  obtained  by  counting  the  branches  of  the  star.  By  means  of  (6)  the 
mass/charge  ratio  can  be  determined  by  adjusting  a  standing  (mostly  eight  branch)  star  pattern 
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and  by  measuring  the  voltage  and  the  frequency  of  the  corresponding  drive. 

Illuminating  the  negatively  charged  particle  with  very  faint  UV  light  the  electrons  will 
be  removed  from  the  particle  one  by  one^.  Some  second  succession  can  be  readily  adjusted  by 
using  faint  enough  UV  light.  If  the  particle  to  be  analyzed  is  in  a  standing  star  pattern  the  loss 
of  electrons  will  cause  the  pattern  to  spin  again.  Then  the  pattern  can  be  stabilized  by  reducing 
the  drive  frequency.  The  mass/charge  ratios  before  and  after  the  electron  loss  can  be  expressed 
by  means  of  (6): 

m  _  n  V  m  ^  ^ 

Ze~  2  ^2  qJ  iZ-l)e  '  2 

where  Z  is  the  ionization  level,  e  is  the  elemental  charge  1.6  x  10'^®  as,  Q, 
drive  frequencies.  The  absolute  mass  can  be  calculated  from  (7): 
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and  Q2  ^re  neighbor 


Since  the  particle  carrier  at  least  several  hundred  unit  charges  the  measurement  can  be  repeated 
a  number  of  times,  thus  making  the  method  more  dependable. 

Laser  Fragmentation/Ionization  MS 

The  LFI  mass  spectrum  of  a  1.0  pirn  polystyrene  particle  illustrated  in  Figure  4  is 
characterized  by  the  occurrence  of  mass  peaks  at  every  nominal  m/z  value  within  the  selected 
mass  range.  Furthermore,  the  mass  intensities  are  quite  high  (tjqjically  10^-10^  "ion  counts  full 
scale).  Finally,  the  spectrum  appears  to  be  dominated  by  organic  compound  signals  in  view  of 
the  frequent  occurrence  of  14,  15  or  28  amu  intervals  rather  than  the  m/z  23  or  39  intervals 
characteristic  of  sodium  and  potassium  ion  contributions  or  the  12  amu  intervals  associated  with 
carbon  cluster  formation.  Unfortunately,  detailed  chemical  interpretation  of  the  mass  spectrum 
in  Figure  3  is  presently  impossible  since  the  signals  do  not  directly  correspond  to  the  expected 
styrene  and  styrene-like  decomposition  products.  The  most  likely  explanation  is  that  the  high 
initial  concentrations  of  ions  in  the  ion  trap  are  causing  significant  shifts  in  the  mass  scale 
calibration.  Similarly,  the  B.  subtilis  spore  spectrum  in  Figure  5  also  provides  high  ion 
intensities  combined  with  a  current  lack  of  biochemical  interpretability.  Complex  secondary 
reactions,  e.g.,  due  to  excessive  laser  intensities  and/or  an  unknown  mass  scale  calibration  shifts 
due  to  space  charge  phenomena  in  the  trap  may  well  be  responsible  for  these  problems.  Before 
any  detailed  conclusions  can  be  drawn,  however,  it  will  be  necessary  to  confirm  the  organic 
nature  and  information  content  of  the  ion  signals  observed.  Since  the  ion  trap  system  can  be 


operated  in  a  tandem  MS  mode  it  should  be  possible  to  produce  daughter  ion  spectra  and  thus 
identify  some  of  the  major  peak  signals  in  Figures  4  and  5. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions: 


1)  Micron-  and  submicron-sized  particles  can  be  trapped  from  electrospray  sources  from 
vaporized  suspensions  in  H2O,  or  directly  from  air. 

2)  Trapped  particles  can  be  pumped  down  to  high  vacuum  (<10'^  torr)  conditions  without 
loss  of  stability. 

3)  Under  high  vacuum  conditions  particles  describe  a  star  like  pattern  in  the  equatorial 
plane  that  can  be  used  for  precise  mass/charge  ratio  determinations. 

4)  Use  of  a  Nd  YAG  laser  beam  at  1.06  [xm  produces  intense  ion  signals  from  single 
styrene  particles  or  B.  subtilis  spores  which  appear  to  represent  complex  organic 
compounds. 

Recommendations: 


1)  Further  experiments  are  needed  to  optimize  Nd  YAG  laser  pulse  intensities  beam 
alignment  and  choice  of  wavelengths  (the  1.06  |xm  wavelength  can  be  doubled,  tripled 
or  quadrupled),  especially  in  order  to  increase  shot-to-shot  reproducibility. 

2)  Use  of  tandem  MS  techniques  should  enable  identification  of  the  presently  unknown  ion 
signals. 

3)  A  systematic  of  different  types  of  aerosol  particles,  supported  by  pattern  recognition 
methods  will  be  needed  to  evaluate  the  practicality  of  the  LFI-MS  technique. 
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Proceeding  from  the  pseudopotential-well  approximation  for  ion  motion  in  a  quadnipole  ion 
trap,  mathematical  expressions  are  derived  to  describe  the  excitation  amplitude  of  an  ion 
packet  at  a  given  mass-to-charge  ratio.  Ion-neutral  collisions  are  incorporated  to  describe  the 
damping  of  ion  trajectories  and  to  describe  the  distribution  of  individual  ion  trajectories 
about  a  mean  amplitude  for  the  ion’ packet  The  rate  of  increase  of  the  amplitude  during 
scanning  is  related  to  expressions  that  describe  the  amplitude  dispersion  of  the  ions  at  the 
time  of  ejection  from  the  trap,  which  is  operating  in  a  resonance  ejection  scanning  mode  to 
describe  the  temporal  line  width  of  the  ejected  ion  packet  The  temporal  line  width  is  related 
to  mass  resolution  under  a  number  of  different  scanning  conditions.  Included  in  the 
discussion  are  considerations  of  the  effect  on  resolution  of  the  resonance  excitation  voltage, 
temperature,  pressure,  noise,  and  buffer-gas  composition.  An  expression  for  the  maximum 
possible  resolution  at  high  ion  mass-to-dharge  ratios  is  developed,  and  these  results  are 
compared  to  an  existing  theoretical  construction.  The  expressions  derived  under  the  pseu- 
dopotential-well  approximation  are  further  extended  to  high  q,  values  and  compared  to 
experimental  data  previously  published  by  two  other  researchers,  (J  Am  Soc  Mass  Spectrom 
1994.  5,  676-688) 


Recent  work  on  high  mass,  high  resolution  ion 
trap  mass  spectrometry  has  shown  that  three- 
dimensional  Paul-type  [1]  quadnipole  instru¬ 
ments  are  capable  of  resolution  exceeding  one  million 
[2,  3]  and  detection  of  ions  with  mass-to-charge  ratios 
greater  than  50,000  [4]  during  operation  in  the  reso¬ 
nance  ejection  mode  [5].  Although  the  theory  that 
describes  the  mass  scanning  properties  of  ion-trap-type 
instruments  in  both  the  mass-selective  instability  mode 
[6,  7]  and  in  the  resonance  ejection  modes  [8]  is  well 
established,  even  in  systems  incorporating  nonlinear 
hexapole  and  octapole  fields  [9],  theory  that  describes 
the  resolution  characteristics  for  the  purpose  of  further 
optimization  of  instrument  performance  has  been  slow 
to  develop. 

In  a  recent  article,  Goeringer  et  al.  [10]  proposed  a 
theoretical  basis  for  understanding  high  resolution 
mass  spectra  obtained  in  quadnipole  ion  traps  via 
resonance  ejection.  This  theory  is  based  on  two  funda¬ 
mental  assumptions.  The  first  assumption  requires  that 
a  pseudopotenfial-weU  description  of  ion  motion, 
which  was  extended  to  include  ion-molecule  colli¬ 
sions,  allows  the  amplitude  excitations  of  the  ion  to  be 
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described  by  the  response  of  a  damped  harmonic  oscil¬ 
lator  to  an  external  driving  frequency.  The  basis  of  this 
assumption  is  the  work  of  Major  and  Dehmelt  [11]. 

By  using  this  assumption,  the  authors  derived  the 
amplitude  response  of  an  ion  in  the  case  of  a  Linearly 
scanned  excitation  frequency.  This  response  showed  a 
clear  dependence  on  scan  rate  and  ion  damping,  and 
the  width  of  this  response  line  in  frequency  space  was 
subsequently  related  to  a  mass  line  width  to  determine 
the  mass  resolution.  The  use  of  this  frequency  line 
width  for  ion  excitations  implies  a  second  assumption. 
This  assumption  is  that  frequency  dispersion  of  the 
amplitude  excitation  is  equivalent  (or  at  least  propor¬ 
tional)  to  the  spatial  dispersion  of  the  ion  packet  at 
ejection.  In  effect,  this  assumption  implies  that  ions  of 
a  given  resonant  frequency  are  being  ejected  from  the 
trap  during  the  entire  time  that  the  excitation  ampli¬ 
tude  of  the  ion  packet  is  nonzero,  so  that  the  temporal 
line  width  is  the  same  as  if  the  ions  were  not  ejected  at 
all. 

In  addition  to  not  directly  addressing  the  physics  of 
ion  packet  ejection,  this  second  assumption  is  trou¬ 
bling  because  it  does  not  allow  for  a  description  of  the 
effect  of  excitation  amplitude  on  mass  resolution.  Fxir- 
thermore,  it  is  well  known  that  a  classical  externally 
driven  damped  harmonic  oscillator  has  no  implied 
amplitude  dispersion  [12]  beyond  that  implied  by  tran¬ 
sient  excitations  and  phase-angle  differences.  If  scanned 
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sufficiently  slowly  and  without  any  other  amplitude 
dispersion  mechanism,  a  damped  harmonic  oscillator 
with  a  fixed  ejection  amplitude  should  be  capable  of 
infinite  resolution 

Figure  1  clarifies  this  point.  The  wide  distribution  in 
Figure  1  describes  the  ion  packet  oscillation  amplitude 
as  a  function  of  time,  assuming  the  behavior  of  a 
damped  harmonic  oscillator.  The  inset  distribution  near 
the  Zq  axis  represents  an  actual  distribution  of  ion 
oscillation  amplitudes  around  the  mean  at  the  moment 
of  ion  packet  ejection.  In  a  first-order  approximation, 
the  actual  temporal  line  width  at  ejection  (inset  near 
the  art  axis)  is  governed  by  the  amplitude  dispersion 
of  the  ion  trajectories  and  the  rate  of  amplitude  in¬ 
crease.  This  line  width  may  be  considerably  thinner 
than  the  temporal  line  width  of  the  amplitude  excita¬ 
tion  of  ions  that  are  not  ejected  from  the  trap.  In  fact,  if 
this  amplitude  dispersion  approaches  zero,  the  tempo¬ 
ral  line  will  also  approach  zero,  corresponding  to  infi¬ 
nite  resolution. 

Furthermore,  this  second  assumption  is  unnecessary 
because  the  same  mathematical  development  of  ion- 
neutral  momentum  transfer  processes  used  to  describe 
the  collisional  damping  of  ion  trajectories  also  can  be 
used  to  describe  the  amplitude  dispersion  of  the  ion 
packet  at  the  point  of  ejection  [13].  The  amplitude 
dispersion  of  the  packet  then  may  be  related  to  the  rate 
of  increase  in  the  ion  packet  oscillation  amplitude  at 
ejection  to  obtain  temporal  and  mass  resolution  line 
widths.  These  line  width  expressions  yield  resolution 
descriptions  that  may  be  compared  to  the  resolution 
expressions  obtained  from  the  previous  frequency  line 
width  assumption  and  to  various  experimental  mea¬ 
surements. 


Theory 

As  illustrated  in  Figure  1,  the  temporal  line  width  of  a 
group  of  ions  ejected  at  a  given  mass-to-charge  ratio 


OTt 


Figure  1.  The  relationship  between  ion  packet  excitation  (full 
width  curve),  amplitude  dispersion  of  the  ion  packet  (inset  left), 
and  the  actual  temporal  line  width  of  the  mass  peak  represented 
by  the  ion  packet  (inset  center). 


may  be  approximated  by  the  expression 

AA 

^  dA/dt 


(1) 


Here  A  f  is  the  temporal  line  width,  A  A  is  the  ampli¬ 
tude  dispersion  at  ejection,  and  dA/dt  is  the  rate  of 
increase  of  the  ion  packet  amplitude  at  ejection.  By 
using  the  pseudopotential- well  approximation  [7,  11], 
we  compute  dA/dt  on  the  basis  of  considerations 
similar  to  those  used  by  Goeringer  et  al.  [10],  whereas 
A  A  is  derived  from  dispersion  considerations  in 
ion-molecule  collisions. 


Axial  Excitation  of  Ion  Motion 

Ignoring  the  radial  motion  equations,  we  describe  the 
axial  motions  as 

d^z  dz  . 

■^  +  ^  =/s  (2) 

where 

q.Ci 

v/2mr, 

ieV 

The  terms  are  defined  as  follows:  z  is  the  mean  axial 
ion  position,  V  is  the  fundamental  rf  drive  voltage,  e  is 
the  ion  charge,  m  is  the  ion  mass,  is  the  ring 
electrode  radius,  O  is  the  fundamental  angular  drive 
frequency,  is  the  supplemental  end  cap  excitation 
voltage,  (o^  is  the  supplemental  end  cap  excitation 
frequency,  and  c  is  the  reduced  collision  frequency 
discussed  below.  The  term  is  the  secular  frequency 
of  unforced  resonant  oscillations,  and  f  is  amplitude 
of  external  forcing  normalized  per  unit  mass. 

The  substitution  of  z  =  Zit)exp(io)J)  into  eq  2 
along  with  an  assumption  that  Z(0  is  a  slowly  vaiying 
function  of  time,  a>  ~  and  co  >  c  allows  us  to  write 


where  =  co^  —  co. 

It  may  be  shown  that  if  the  secular  frequency  is 
swept  linearly  with  time  such  that  A  ca  =  at,  where  a  is 
the  scan  rate  in  radians  per  square  second,  then  Z  can 
be  written  as 

2  (ijfa  ^  T  2  j 


678  ARNOLD  ET  AL. 


J  Am  Soc  Mass  Spectrom  1994,  5,  676-688 


where  x  =  and  r  =  2/c,  The  amplitude  of  the 
secular  oscillation  is  determined  by  the  absolute  value 
of  Z(:rX  such  that  the  amplitude  A  ^  |2(x)|.  To  sim¬ 
plify  subsequent  operations,  we  define  F(x)  as  a  di¬ 
mensionless  amplitude 


2  ojy/a  2  cj/a 

Fix)  =  — _|Z(:c)|  =  ——Aix)  (5) 

/s  /i 


where  Fix)  is  given  by 


Direct  evaluation  of  these  •  integrals  is  difficult,  and 
closed-form  solutions  do  not  exist,  but  as  shown  in  refs 
10  and  11,  it  is  possible  to  utilize  approximations  for 
the  two  cases  where  >  1  and  <  1.  These 
integrals  also  may  be  evaluated  numerically  in  regions 
where  these  simple  approximations  are  not  available. 

For  <  1,  we  utilize  a  result  obtained  by  inte¬ 
gration  by  parts,  so  that 


Fix) 


Vl  +  ar^x^ 


(7) 


gives  the  excitation  amplitude.  This  expression  may  be 
differentiated  with  respect  to  time  to  obtain 

dF  i4ar)^  X 

■—  =  -  _  3/2  (S) 

This  latter  equation  describes  the  rate  of  increase  of  the 
ion  envelope  near  ejection. 

For  the  case  of  >  1  we  require  an  approxima¬ 
tion  that  describes  the  rising  side  of  the  amplitude 
envelope  so  that  the  slope  in  the  vicinity  of  resonant 
frequency  at  x  =  0  may  be  determined.  In  this  limit, 
we  can  rewrite  eq  6  in  terms  of  Fresnel  integrals  [14], 
so  that 


shows  a  square  root  dependence  of  the  diffraction 
intensity  pattern  produced  by  an  infinite  straight  edge 
in  wave  optics  [15].  Subsequent  differentiation  with 
respect  to  x  gives 

dF 

dx 


(10) 


Amplitude  Dispersion 

The  average  rate  of  momentum  loss  due  to  collisions  is 
related  to  the  reduced  collision  frequency  c.  A  mathe¬ 
matical  derivation  for  c  is  given  by 


c  = 


M 

m  -r  M 


(11) 


where  m  is  the  ion  mass,  M  is  the  mass  of  colliding 
neutrals,  n  is  the  number  density  of  neutrals,  v^.  is  the 
average  relative  speed  of  the  colliding  ion  and  neutral, 
and  cTp  is  the  collision  cross  section,  and  proceeds 
from  momentum  transfer  arguments  that  are  outlined 
elsewhere  [13],  wherein  the  ion  is  assumed  to  be 
"drifting"  in  a  constant  electric  field.  The  present  de¬ 
velopments  assume  that  although  the  ions  undergo 
accelerated  motion  in  the  ion  trap  oscillations,  the 
average  momentum  damping  is  still  the  same. 

The  same  arguments  that  allow  the  development  of 
the  reduced  collision  frequency  expressioris  also  may 
be  applied  to  velocity  dispersion  of  the  ion  packet  [13] 
caused  by  collisions  with  neutrals.  This  dispersion  is 
based  on  the  discontinuous  nature  of  the  ion  damping. 
Following  Mason  and  McDaniel  [13]  and  considering 
only  the  axial  excitations  of  the  ions,  we  write  the 
relative  velocity  dispersion  as  for  ions  traveling  at  drift 
velocity 


iAV,f 


C*  4m  +  M  M 
3  2m  -h  M  m 


(12) 


where  ( A  7^)^  is  the  velocity  dispersion  about  the  mean 

in  the  z  direction,  is  the  thermal  average  of  the 
square  of  ion  velocity  in  the  z  direction,  k  is  the 
Boltzmann  constant,  T  is  absolute  temperature,  and 
C*  is  a  constant  that  has  a  value  on  the  order  of  1 
depending  on  the  m/M  ratio  and  the  nature  of  ion- 
neutral  scattering. 

As  mentioned  before,  eqs  11  and  12  are  both  devel¬ 
oped  for  the  case  of  constant  electric  field.  Use  of  eq  11 
in  the  present  case  requires  the  assumption  that  the 
overall  ion  packet  is  in  dynamic  equilibrium  with  the 
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neutrals  via  collisions,  so  that  c  may  represent  the 
average  damping  of  the  ion  packet.  Equation  12  repre¬ 
sents  the  fluctuations  around  this  equilibrium,  and  the 
earlier  assumption  that  o)  >  c  implies  that  the  velocity 
dispersion  may  be  averaged  over  time  by  using  the 
root  mean  square  (rms)  value  of  v^.  Further,  because 
that  particle  is  undergoing  simple  harmonic  motion, 
the  expression  v^  —  coz  allows  us  to  relate  the  velocity 
dispersion  to  the  amplitude  dispersion  of  the  ion 
packet. 

The  amplitude  dispersion  of  the  ion  packet  in  dy¬ 
namical  equilibrium  with  the  surrounding  neutrals 
may  therefore  be  represented  as 


1 

AA  -  — 


C*  4m  +  M  M  ^ 
T  2m  +  M 


(13) 


where  is  the  rms  velocity  amplitude  due  to  the 
external  driving  field.  At  the  point  of  ejection,  A^j  = 

substitution  into  eq  13  gives 


AA 


C*  4m  -h  M  M  2 
12  2m  H-  M  In 


(14) 


.  This  expression  indicates  that  the  amplitude  disper¬ 
sion  at  ion  ejection  is  directly  related  to  the  effective 
ioiii  temperature  in  the  axial  direction  [13].  Additional 
consideration  of  velocity  dispersion  in  the  radial  direc¬ 
tion  would  allow  m(Ay)^  (similar  to  eq  12)  to  be 
related  to  the  effective  ion  temperature  as  discussed  in 
Mason  and  McDaniel  [13].  This  fact  has  additional 
implications  for  instrument  resolution  and  spectral 
quality,  which  may  be  degraded  by  collision-induced 
dissociation  of  ions  caused  by  scanning  of  the  ions 
from  the  trap. 


Minimum  Ejection  Line  Width 

Direct  substitution  into  eq  1  of  eq  14  along  with  eq  8  or 
10  yields  an  expression  for  temporal  line  width  At 
depending  on  the  actual  time  of  ion  packet  ejection 
Although  this  substitution  may  be  performed  at  any 
value,  tiie  primary  interest  is  in  the  maximum  obtain¬ 
able  resolution  at  a  given  scan  speed  that  will  occur  at 
the  minimum  of  At  or  the  maximum  of  dA/dt. 

At  the  point  in  time  (=  where  the  ion 

package  is  ejected,  the  amplitude  of  oscillation  is  given 
by  A(Xg|)  =  7^/2^/“.  Based  on  eq  5,  this  ejection  point 
is  determined  by  the  excitation  amplitude  /g,  so  that 


F(x,p 


(15) 


describes  the  excitation  required  to  eject  the  ion  packet 
at  time  x^. 

To  evaluate  the  maximum  of  dA/dt,  eq  5  is  differ¬ 
entiated  with  respect  to  t  whereas  eq  15  is  substituted 


for  /s  to  obtain 

dA  rja  dT{x^)  rja  dlnF(x„) 

-  =  - ; -  =  — •}=■ - ; -  Vlo; 

dt  v^F(Xej)  dx  n  dx 


dA/dt  is  proportional  to  the  ratio  of  dE/dx  to  F.  This 
ratio,  also  expressed  as  d  In  E /dx,  is  plotted  for  vari¬ 
ous  values  of  in  Figure  2.  Two  different  scalings 
are  provided  in  this  figure  because  the  limiting  case 
<<  1  has  a  natural  scaling  of  art,  as  seen  in  eqs  7 
and  8,  whereas  the  case  1  scales  to  as 

seen  in  eqs  9  and  10.  At  a  given  scan  speed  and 
damping,  the  maximum  resolution  depends  on  maxi¬ 
mizing  d]nE/dx,  This  maximum  occurs  at  t^=^ 
- 1  /ar,  as  illustrated  in  Figure  2a,  for  all  scan-rate  and 
damping  values!  This  result  establishes  the  temporal 
shift  of  the  ion  ejection  point  (obtained  at  maximum 
resolution)  relative  to  the  resonant  frequency  of  ion 
oscillations. 

The  maximum  of  dA/dt  at  ion  ejection  as  a  function 
of  and  r  may  be  obtained  for  the  two  limiting  cases 
on  the  basis  of  a  similar  evaluation  as  was  performed 
to  obtain  eqs  8  and  10.  For  a^^^r  <  1,  the  rate  of 


ait 


Figure  2.  Derivative  of  the  log  of  f  ( x)  for  various  scan-rate  and 
damping  factors  as  a  function  of  art  and  All  curves  share 
a  maximum  at  art  =-  -1.  The  maximum  on  each  curve  indicates 
the  optimum  excitation  amplitude  and  resolution. 
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Mass  Dependence  of  Instrument  Resolution 

Equations  24  and  25  describe  the  mass  dependence  of 
the  maximum  obtainable  resolution,  but  in  fact  both  a 
and  T  have  mass  dependence  that  should  be  made 
explicit.  This  affects  both  the  resolution  expressions  of 
eqs  22  through  25  as  well  as  appropriate  axial  excita¬ 
tion  amplitudes  discussed  subsequently  in  eqs  37 
through  40. 

The  mass  dependence  of  a  is  straightforward:  a  is 
directly  related  to  the  mass  scan  rate  via  the  differen¬ 
tial  form  of  eq  21  as  seen  in 


Like  r  also  has  mass  dependence  because  r  =  2/c 
and  c  is  defined  in  eq  11.  In  addition  to  the  explicit 
terms,  and  cr^  may  also  carry  some  mass  depen¬ 
dence.  The  special  case  utilized  by  Goeringer  et  al.  [10] 
involved  an  induced  dipole-interaction  model  between 
the  ion  and  the  neutral  molecules,  because  for  this 
model,  the  collision  frequency  is  truly  independent  of 
ion  velocity  in  In  such  a  case  (by  using  m  >  M)  we 
obtain 


d(o  eoj  d{m/e) 

dt  m  dt 


(26) 


Although  the  frequency  scan  rate  for  a  given  mass-to- 
charge  ratio  remains  constant  over  the  course  of  the 
scan,  the  effective  (frequency)  scan  rate  decreases  with 
increasing  mass.  This  important  result  may  be  coupled 
immediately  to  eq  23  (high  scan  speed  limit)  to  obtain 


v'ars.l  (27) 


This  result,  where  resolution  is  directly  proportional  to 
mass,  may  be  used  to  illustrate  clearly  why  unit  mass 
resolution  is  readily  maintained  over  the  entire  scan 
range  of  the  "classic"  ion  trap  (16,  17]  system. 

A  brief  qualitative  argument  illustrates  this  fact.  It 
is  well  recognized  that  an  ion  trap  that  uses  axial 
modulation  [8]  qualitatively  gives  resolution  behavior 
similar  to  the  classic  mass-selective  instability  scan, 
whereas  various  researchers  [10,  18,  19]  have  noted 
that  line  width  behavior  at  higher  values  (i.e., 
q^  >  0.4)  should  yield  qualitatively  similar  characteris¬ 
tics.  At  the  high  q^  values  (near  q^  =  0.9)  that  are  used 
for  axial  modulation,  JcT  <  C*Mr^(o^/6,  so  that  in  clas¬ 
sic  mass-selective  instability  at  high  scan  speeds,  eq  28 
pertains: 


(28) 


The  ion-mass  m  and  neutral-molecule  buffer-mass 
M  [20]  dependence  of  this  expression  are  in  good 
agreement  with  the  known  behavior  of  the  classic 
mass-selective  instability  scan.  The  scan-rate  depen¬ 
dence  also  agrees  with  high  scan  speed  behavior  de¬ 
scribed  by  Louris  et  al.  [19]  (see  Figure  6a).  This  result 
further  implies  that  resolution  is  independent  of  pres¬ 
sure  if  ions  are  sufficiently  cooled  prior  to  spectrum 
scanning  [6,  7,  21]. 


where  is  the  polarizability  of  the  neutral  and  is 
the  permittivity  of  free  space. 

The  velocity  independence  of  c  is  essential  to  the 
solution  of  eq  2  as  presented  in  eqs  7  through  10.  In 
this  ion-induced  dipole-interaction  model,  this  velocity 
independence  is  preserved,  although  as  seen  in  eq  11, 
the  reduced  collision  frequency  c  in  general  depends 
on  the  relative  velocity  between  the  ion  and  a  colliding 
neutral  molecule.  However,  even  in  a  more  general 
ion-neutral  interaction  model,  the  relative  velocity  Vj, 
(and  therefore  c)  is  still  velocity  independent  if  kT  > 
C*Mr-a>“/6.  This  condition  is  the  same  as  that  which 
will  be  discussed  subsequently  as  a  temperature- 
limited  resolution  regime  (see  eq  32). 

Nevertheless,  this  expression  has  some  general  ap¬ 
plicability  for  lower  mass  ions,  and  utilizing  eq  29,  one 
may  rewrite  eq  22  as 


(30) 


This  is  the  low  scan  speed  analog  of  eq  27. 

It  must  be  noted  that  eqs  27,  28,  and  30  all  show  a 
difference  in  resolution  between  ions  of  the  same 
mass-to-charge  ratio  when  these  ions  have  different 
charge  states.  More  explicitly,  each  equation  can  be 
written  as  a  function  of  m/e  multiplied  by  In 

this  situation,  the  improvement  in  resolution  for  a 
multiply  charged  ion  at  a  given  m/e  value  is  propor¬ 
tional  to  the  number  of  elemental  charges.  It  is  impor¬ 
tant  to  note  that  this  implies  that  under  a  fixed  set  of 
scanning  conditions,  the  maximum  mass  resolution 
will  be  obtained  for  the  singly  charged  ion. 

The  primary  means  of  delineating  the  two  operating 
regimes  has  been  to  consider  the  value  of  relative 
to  1.  This  value  also  has  an  explicit  mass  dependence, 
based  on  eqs  26  and  29,  as  is  given  by 


It  is  clear  that  the  primary  dependence  in  this  expres¬ 
sion  is  pressure,  whereas  all  other  factors  are  ex- 
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pressed  as  the  square  root.  Nevertheless,  under  full 
scan  experiments,  the  mass  dependence  does  provide 
for  differences  of  a  factor  of  3  or  4  for  different  ions  in 
a  given  scan,  yielding  distinctive  resolution  behavior. 


Temperature  and  Collisional  Focusing 

The  effect  of  temperature  on  resolution  arises  from  the 
ion  packet  dispersion  term  in  the  denominator  of  eqs 
22  through  25.  The  dependence  of  resolution  on  tem¬ 
perature  is  governed  by  the  relative  values  of  the 
temperature  itself  and  the  frequency  of  supplemental 
excitation.  The  critical  point  occurs  when 


kT^ 


T 


Mr^co^ 


(32) 


where  is  the  mean  thermal  velocity  of  the  neutrals. 
This  point  occurs  when  the  average  axial  ion  velocity 
at  ejection  {v^  =  exceeds  the  mean  thermal 

velocity  of  the  neutrals.  For  ejection  frequencies 
values)  below  this  critical  point,  the  resolution  is  lim¬ 
ited  by  the  dispersion  of  the  ''stationary"  ion  packet 
before  excitation,  which  is  a  function  of  system  tem¬ 
perature  and  may  be  represented  as 


m  1 

Am  ^ 


2 


(33) 


This  expression  clarifies  that  increases  in  ejection  fre¬ 
quency  greatly  enhance  the  resolution  until  the  condi¬ 
tion  of  eq  32  is  met.  In  a  standard  commercial  ion  trap 
mass  spectrometry  (ITMS)  instrument  (r^  =  1  cm,  H  == 
6.911  X  10^  s“^)  at  0°C,  the  condition  of  eq  32  is 
computed  to  occur  at  a  supplemental  frequency  of  29.3 
kHz  =  0.075)  for  helium  buffer  gas  and  11.1  kHz 
(q^  -  0.028)  for  nitrogen  (assuming  C*  =  1).  Above 
these  values,  the  increase  in  resolution  with  increasing 
frequency  will  be  a  modified  linear  relationship  de¬ 
scribed  by 


m  (w/i/8)V'6m7cW 
Am  Tj-ira  /2  +  1/r 


C* 

6 


(34) 


This  is,  of  course,  the  opposite  linuting  case  from  eq 
33.  In  this  case,  we  see  that  buffer  molecular  mass 
becomes  an  important  limiting  factor  in  higher  ejection 
frequency  cases.  In  either  case,  increasing  the  ion  mass 
or  the  excitation  frequency  yields  increased  resolution. 

Figure  4  illustrates  the  relationship  between  tem¬ 
perature,  buffer  mass,  and  excitation  frequency  for  five 
different  buffer  gases  in  both  low  and  high  scan  speed 
limits.  At  low  temperatures  (the  buffer-mass-limited 
case  of  eq.  28),  the  condition  of  maximum  resolution  is 
obtained,  whereas  increasing  temperature  provides  a 
degraded  resolution  although  the  inflection  point  may 


Log[kT/(r,co)2] 


Log[kT/(r^o>)2] 

Figure  4.  Temperature  dependence  of  resolution  for  various 
buffer  gases:  (a)  high  scan  speed  limit  r  >  1;  (b)  low  scan 
speed  limit  r  ^  1. 


be  varied  by  adjustment  of  the  excitation  frequency. 
One  interesting  feature  of  the  low  scan  speed  limit  is 
the  fact  that  helium  rather  than  hydrogen  provides  the 
highest  resolution.  This  is  primarily  due  to  the  signifi¬ 
cantly  lower  polarizabilityfby  more  than  a  factor  of  2) 
of  helium. 

The  temperature  of  neutral  buffer  molecules  has 
been  used  in  the  development  of  amplitude  dispersion 
equations  12  through  14.  This  approach  is  consistent 
with  the  pseudopotential-well  model  developments  of 
Major  and  Dehmelt  [11].  It  is  clear  that  the  actual  ion 
temperatures  reported  in  the  literature  [21-23],  which 
are  typically  higher  than  the  buffer  ion  temperature, 
depend  also  on  the  effective  ion  temperatures  pro¬ 
duced  by  non-thermal  excitations  of  the  ion  trajecto¬ 
ries.  Several  means  for  these  excitations  have  been 
discussed  including  nonlinear  resonances  [9,  24],  drive 
frequency  harmonics  [24],  parasitic  electrical  signals, 
me^anical  vibrations,  and  electrical  noise  [25]. 

Although  these  excitation  sources  must  be  recog¬ 
nized  for  their  temperature-related  effects  in  any  ion 
chemistry  experiment,  it  is  clear  that  the  hottest  point 
for  ions  in  virtually  all  ion  trap  experiments  occurs  at 
ion  ejection.  At  ejection,  the  amplitude  dispersion  de¬ 
scribed  by  eq  14  has  a  direct  thermal  component  and 
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an  axial  excitation  component.  The  axial  excitation 
component  of  this  dispersion  is  proportional  to  the 
product  A(o,  as  was  described  in  the  derivation  of  eq 
14,  for  each  frequency  component.  As  long  as  the  Aoj 
product  for  other  non-secular  oscillation  frequencies  is 
small  with  respect  to  the  dispersion  factors 

for  these  other  oscillations  may  be  neglected  in  the 
spatial  amplitude  dispersion  term. 


Excitation  Amplitude  and  Mass  Shift 

During  the  development  of  minimum  ejection  line 
width  expressions,  eq  15  described  the  relationship 
between  the  ejection  point  of  the  ion  packet  and  the 
supplementary  excitation  amplitude.  Further  it  was 
observed  that  dA/dt  has  a  clear  maximum  for  given  a 
and  T  values  when  the  ejection  point  occurs  at  At^j^  as 
given  by 


(35) 

(see  Figure  2).  This  time  is  measured  relative  to  the 
resonance  point  in  the  scan  (e.g.,  t  —  0  occurs  at  a>  = 
Wg).  Substituting  eqs  26  and  29  into  this  expression 
gives  eq  36,  which  illustrates  the  effect  of  pressure, 
scan  rate  {d{m/e)/dt),  and  ejection  frequency  on  the 
ejection  time  at  maximum  resolution: 

-1-/1-^))/?  - 

This  ejection  time  implies  a  mass  shift  that  can  be 
determined  based  on  eq  21  to  offer  a  comparison 
between  mass  shift  and  mass  resolution.  The  relative 
mass  shift  at  maximum  resolution  is  given  by 

m 

Am 

which  is  independent  of  scan  rate,  but  clearly  depen¬ 
dent  on  ion  mass-to-charge  ratio,  buffer  mass,  and 
pressure.  For  low  scan  speeds  (a^^V  <  1)  this  value 
may  be  directly  compared  to  those  in  eqs  22  and  30. 

This  result  is  somewhat  limited  in  that  it  only 
applies  to  the  mass  shift  after  the  excitation  amplitude 
/s  is  adjusted  to  obtain  maximum  resolution.  If  the 
amplitude  of  the  secular  drive  is  not  adjusted  to  yield 
maximum  resolution  for  the  ion  mass  of  interest,  it  is 
clear  that  eq  15  will  govern  the  mass  shifts.  The  maxi¬ 
mum  delay  possible  in  an  ion  ejection  system  occurs  at 
the  maximum  value  of  amplitude  envelope.  In  Figure 
2  this  corresponds  to  the  point  where  d  In  ¥/dx  =  0. 
At  scan-rate  extremes,  ejection  can  only  be  delayed  for 
a  small  time  after  resonance.  For  <  1,  the  maxi¬ 
mum  delay  occurs  approximately  at  t  ==  -f-2T,  whereas 
for  >  1,  the  maximum  delay  occurs  at  about 
t  =  +2/a^^^.  By  overdriving  the  axial  excitation  am- 


(jjm 


OJT  = 


irne  y 


(37) 


plitude,  ions  may  be  ejected  at  any  time  prior  to 
resonance,  so  that  the  minimum  ejection  time  t  is  not 
Umited  by  any  boundary. 

Because  maximum  resolution  occurs  at  — 
as  seen  in  Figure  2,  substitution  of  this 
value  into  eq  15  provides  an  expression  for  optimum 
excitation  amplitude.  Combined  with  the  defiiiition  of 
/s  in  eq  2,  f  can  be  directly  related  to  the  excitation 
voltage  amplitude,  so  that 


Icomr^fa 

^  H-1/^’t) 


(38) 


The  behavior  of  as  a  function  of  is  the  inverse 
of  the  behavior  of  resolution  m  Figure  3.  The  value  of 
Vg  may  be  approximated  in  a  manner  similar  to  that 
for  m/Am  in  eqs  24  and  25.  The  limiting  cases  are 
given  by  the  expressions 


2i/2 


comr^ 


a 

'tr 


and 


iflcomr^  /  1  \ 

e  [  r  j 


(39) 


(40) 


These  expressions  may  again  be  substituted  with  eqs 
26  and  29  to  illustrate  the  mass  dependencies  at  maxi¬ 
mum  resolution  conditions.  Further,  substituting  eq  26 
for  the  case  >  1  yields 


V:  ==  2fI<or^ 


mcj  dim/e) 


Tre  dt 

whereas  substituting  eq  29  in  the  case 


vCt>1  (41) 
■*  1  yields 


Vg  =  ifl  (O'Trnr^^  ^  ■  far  1  (42) 

Of  particular  interest  is  the  mass  dependence  of  both 
expressions,  which  indicate  that  excitation  frequency 
amplitudes  may  benefit  from  scarming  in  a  manner 
similar  to  the  rf  scarming  of  the  instrument.  It  should 
be  noted  that  eq  41  describes  a  behavior  that  is  qualita¬ 
tively  similar  to  that  of  a  classical  ITMS  system  for 
adjustment  of  the  excitation  amplitude  relative  to  the 
scan  rate  and  ion  mass  variables  [26]  and  agrees  weU 
with  ref  19  at  4.2  V. 

This  analysis  is  based  on  the  assumption  that  a 
simple  dipole  field  describes  the  electric  field  pro¬ 
duced  by  the  supplemental  end  cap  voltages.  Actual 
field  values  vary  considerably  from  this  dipole  ideal. 
Even  near  the  origin,  correction  values  for  the  axial 
excitation  electric  field  have  been  reported  as  0.8  [27] 
and  0.877  [28]  based  on  experimental  and  theoretical 
considerations,  respectively.  These  dipole  field  correc- 
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tion  values  are  a  function  of  both  radial  and  axial 
positions,  and  they  are  expected  to  be  considerably 
smaller  with  radial  displacements  near  the  ring  and  to 
increase  to  unity  near  the  center  of  the  end  caps. 

Limitations  and  Comparison  to  Existing  Theory 

The  present  theoretical  construction  is  an  important 
departure  from  previous  theoretical  work  [10]  for  three 
reasons.  First,  the  high  scan  speed  behavior  is  not 
dependent  on  pressure  and  proceeds  as  rather 
than  ar  as  in  the  previous  analysis.  It  should  be  noted 
that  this  result  is  consistent  with  a  large  body  of  data 
that  show  resolution  to  be  constant  over  a  broad  range 
of  pressures  when  sufficient  time  is  allowed  to  ther- 
malize  the  ions  before  scanning  [7,  22]. 

Second,  the  present  theory  predicts,  via  eqs  39 
through  42,  the  excitation  voltages  required  to  obtain 
maximum  resolution  at  any  scan  speed.  However,  we 
reiterate  that  important  nonlinear  field  effects  were 
disregarded  in  the  derivation  of  these  two  equations.  A 
correction  factor  is  necessary  to  predict  required  exci¬ 
tation  voltages  accurately  for  a  given  resolution;  one 
correction  factor  is  probably  obtained  best  from  experi¬ 
ment,  but  the  scan-rate  dependencies  and  temporal 
line  shifts  should  not  be  altered  by  these  considera¬ 
tions. 

Third,  the  ion  dispersion  description  implies  that 
resolutions  are  directly  dependent  on  the  relative 
masses  of  the  buffer  neutrals  and  the  ions  themselves 
owing  to  collisional  focusing.  Furthermore,  this  effect 
is  independent  of  the  scan-rate  regime.  This  effect 
describes  how  higher  resolutions  are  produced  for 
higher  masses  in  any  continuous  scan.  Tliis  term  also 
incorporates  the  effects  of  temperature  (and,  qualita¬ 
tively,  electrical  noise)  into  the  consideration  of  resolu¬ 
tion. 

Figure  5  illustrates  a  predicted  resolution  compari¬ 
son  between  the  present  result  and  the  previous  fre¬ 
quency  line  width  mode  [10].  The  resolution  predicted 
in  the  present  theory  is  actually  higher  in  both  scan-rate 
regimes.  In  the  slow  scan  regime  the  resolution  is 
higher  owing  to  the  collisional  focusing  of  the  higher 
mass  ions  by  the  light  neutrals.  This  focusing  effect  is 
also  present  at  high  scan  speeds,  but  it  is  more  impor¬ 
tant  to  note  that  the  behavior  of  the  resolution 

implies  higher  resolution  even  without  this  focusing 
effect. 

Note  that  several  important  sources  of  amplitude 
dispersion  have  been  neglected  in  the  present  theory  to 
focus  on  the  ion  behavior  in  a  damped  harmonic 
osdliator  model.  During  the  previous  discussion  of 
temperature,  electrical  noise  was  dismissed  as  a  major 
source  of  spatial  ion  packet  dispersion,  but  temporal 
and  frequency  dispersion  effects,  which  could  result 
from  noise  in  the  radiofrequency  drive  or  other  electri¬ 
cal  and  mechanical  sources,  were  not  considered.  It  is 
worth  noting  that  A  is  directly  proportional  to  AV 
in  the  drive  voltage  and  that  this  may  yield  significant 


LogCa^'^T) 


Figure  5.  (a)  Comparison  of  mass  resolution  at  four  different 
m/M  values  in  the  low  temperature  limit.  Both  upper  and  lower 
bounds  are  shown.  The  behavior  of  the  previous  frequency  line 
with  model  is  also  shown  for  comparison,  (b)  Specific  compari¬ 
son  of  the  present  model  at  m/M  =  125  with  the  existing  fre¬ 
quency  line  width  model.  The  upper  four  broken  curves  illus¬ 
trate  the  present  model  at  three  damping  factors,  whereas  the 
lower  3  illustrate  the  frequency  line  width  model  of  the  same 
three  values.  The  two  solid  lines  indicate  the  maximum  obtain¬ 
able  resolution  as  a  function  of  scan  rate  for  the  two  models. 


temporal  fluctuations  in  the  amplitude  of  the  ion  packet 
during  scanning. 

Geometrical  considerations  also  have  been  ne¬ 
glected.  In  a  standard  commercial  instrument,  the  axial 
position  that  corresponds  to  the  ion  exit  aperture  is  a 
function  of  radius.  If  the  exit  aperture  is  only  3.5  nun 
in  radius,  the  A  A  of  the  aperture  may  be  as  large  as 
0.04  r^,  which  is  the  same  as  collision-based  dispersion 
produced  by  an  m/M  ratio  of  92  (e.g.,  m/z  370  in 
He  buffer)  in  the  low  temperature  limit.  Clearly  for 
higher  masses  the  aperture  size  could  easily  be  the 
resolution-limiting  element. 

Geometrical  variations  already  mentioned  with  re¬ 
spect  to  the  supplemental  axial  excitation  field  will 
also  reduce  resolution  because  the  effective  excitation 
voltage  is  a  function  of  radius  in  the  trap.  This  second 
effect  has  the  potential  to  be  more  important  than  the 
first,  given  the  strong  radial  dependence  of  the  excita¬ 
tion  field  strength  already  discussed.  It  is  our  opinion 
that  these  geometrical  effects  are  the  dominant  resolu- 
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tion-liirdting  effects  for  high  mass  analyses  in  commer¬ 
cial  instruments.  A  full  accounting  of  these  geometrical 
effects  is  beyond  the  scope  of  this  article,  but  it  is  clear 
that  both  effects  may  be  mitigated  by  use  of  small 
apertures,  which  limit  ion  detection  to  those  trajecto¬ 
ries  near  the  axis  of  symmetry. 


Implications  for  High  Mass  Analyses 

The  driving  force  behind  examinations  of  ion  trap 
resolution  in  mass-selective  instability  scans  has  been 
attempts  to  obtain  high  resolution  mass  spectra  at  very 
high  mass-to-charge  ratios  [10,  29].  Because  a  given 
instrument  is  generally  limited  to  a  maximum  operat¬ 
ing  voltage  either  by  electrical  discharge  or  engineer¬ 
ing  constraints,  three  means  have  been  used  to  extend 
the  mass  range  of  the  ion  trap  to  allow  ejection  of  high 
mass  ions:  reduction  of  the  resonance  ejection  operat¬ 
ing  point  [5,  30],  reduction  of  the  operating  fre¬ 
quency  n  [31,  32],  and  reduction  of  the  trap  radius 
[31,  32].  When  one  or  more  of  these  means  is  coupled 
with  slow  scan  speeds,  a  combination  of  high  mass 
analysis  and  high  resolution  may  be  obtained. 

Use  of  any  one  of  these  three  means  to  expand  the 
mass  range  results  in  temperature-limited  resolution 
because  kT  >  In  the  case  of  this  limit,  eq 

32  may  be  rewritten  as 


m 

Am 


{slr^{q,af 

\6'iTnefla^MkT 


C* 

far  <l,kT  ^ 

6 

(43) 


which  is  the  primary  resolution  expression  for  high 
mass  analyses. 

The  m^/^  dependence  in  eq  43  indicates  promising 
resolution  possibilities,  but  is  in  itself  misleading.  If 
the  operating  voltage  could  increase  without  limit,  this 
would  describe  the  mass  dependence  of  resolution.  In 
practice,  the  highest  mass-to-charge  ratio  detectable  is 
directly  proportional  to  the  maximum  operating  volt¬ 
age  according  to  the  definition  of  q^  in  eq  2.  Substitu¬ 
tion  of  the  q^  expression  into  eq  43  allows  expected 
resolution  at  high  masses  to  be  written  as 

m 

Am  rrnrlO}yj2(x^mMkT 
C* 

far  <l,kT  >  —MrloA  (44) 

6 

Equation  44  clarifies  that  once  the  maximum  value  of 
V  has  been  reached,  resolution  as  expected  decreases 
with  increasing  mass  because  the  operating  point  q^ 
must  be  reduced  to  extend  the  mass  range  of  the 
instrument. 

It  is  not  clear  that  reductions  in  or  fl  in  eq  44  are 
accompanied  by  an  increase  in  the  operating  point  and 


that  this  procedure  may  only  continue  until  a  maxi¬ 
mum  effective  operating  point  for  the  trap  is  reached. 
This  maximum  effective  value  from  other  experiments 
appears  to  be  q^  -  0.73  [10,  18]  rather  than  the  ex¬ 
pected  stability  limit  of  q^  =  0.91. 

A  clear  impression  of  the  maximum  possible  resolu¬ 
tion  at  high  mass  may  be  obtained  by  a  different 
substitution  of  q^  into  eq  43.  Substituting  for  the  prod¬ 
uct  we  obtain 


Am  47m^2apMkT 


far  <t  l,kT  >  — MrlaA 
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(45) 


which  describes  the  maximum  possible  resolution  for 
an  ion  of  mass  m  at  operating  point  q^.  The  ratio  V /r^ 
comes  from  electrical-discharge  considerations  and 
represents  a  maximum  value  that  is  related  to  a  maxi¬ 
mum  electric  field  strength.  It  is  also  important  to 
notice  that  this  result  is  independent  of  the  electrical 
charge  of  the  ion. 

It  is  important  to  recognize  that  this  result  is  based 
on  the  temperature-limited  regime,  but  that  reductions 
in  operating  frequency  actually  increase  the  secular 
frequency  o)  so  that  the  kT  >  C'*Mrl(D^/6  condition 
may  be  altered.  It  is  also  important  to  note  that  a 
similar  evaluation  of  eq  33  to  describe  the  maximum 
possible  resolution  for  any  ion  is  more  difficult  and 
beyond  the  scope  of  this  paper. 

The  analysis  that  produced  eqs  43  through  45  uti¬ 
lized  the  ion-induced  dipole-interaction  model.  Once 
an  ion  is  sufficiently  large,  the  probability  that  a  neu¬ 
tral  molecule  will  interact  with  the  charge  center  rather 
than  a  neutral  portion  of  the  ion  structure  becomes 
increasingly  small.  The  ion-neutral  interaction  radius  is 
given  by 


b 


2  _ 
o 


m  >  M 


(46) 


When  the  ion  dimensions  are  on  the  order  of  we 
may  expect  that  a  hard-space-type  scattering  model 
will  take  over,  so  that  the  mass  dependence  of  c,  as  per 
eq  11,  will  relate  to  the  molecular  cross  section  of  the 
ion  (i.e.,  (Tq  is  proportional  to  Once  this  limit  is 

reached,  the  analog  on  eq  43  becomes 


A/n  nf^ 


C* 

far  <  \,kT  >  —Mr’^oA 
6 

(47) 


for  very  large  m.  Similar  analogies  to  eqs  44  and  45 
indicate  that  resolution  decreases  with  increasing  mass 
(once  the  voltage  limit  is  reached)  by  factors  of 
and  respectively. 

Ultimately  the  improvement  of  resolution  will  cease 
to  keep  up  with  the  increasing  mass.  By  using  values 
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for  the  polarizability  of  N2  and  He  included  in  Mason 
[13],  we  can  show^that  the  collision  radius  is  12.5  A 
for  N2  and  7.3  A  for  He  with  respect  to  a  singly 
charged  ion.  These  values  are  easily  exceeded  by  the 
dimensions  of  a  singly  charged  ion  at  10,000  u,  so  that 
expected  resolutions  for  such  high  mass  ions  will  re¬ 
quire  extremely  low  pressures  to  compensate  for  the 
less  favorable  hard-sphere-type  interaction  conditions. 

Beyond  the  Psendopotential-Well  Approximation 

At  several  times  during  this  discussion,  we  alluded  to 
behavior  of  resonance  ejection  scans  for  values 
beyond  the  0.4  limit  of  the  pseudopotential-well  ap¬ 
proximation.  This  is  of  course  because  the  resonance 
phenomena  themselves  are  still  active  at  these  higher 
q.  values  [8,  9].  Furthermore,  we  believe  there  is  noth¬ 
ing  conceptually  flawed  about  the  use  of  this  descrip¬ 
tion  at  higher  q^  values  if  appropriate  considerations 
are  developed.  It  is  for  this  reason  that  we  have  typi¬ 
cally  utilized  (o  rather  than  q^Cl  in  resolution  (and 
other)  expressions  so  that  they  may  be  extended  to 
high  q^  values. 

The  first  consideration  that  must  be  addressed  to 
extend  the  present  work  beyond  the  pseudopotential- 
well  region  is  the  use  of  a  harmonic  oscillation  model 
to  determine  the  amplitude  dispersion  in  eq  14.  Out¬ 
side  the  pseudopotential  region,  17^5  as  used  in  eq  13 
is  no  longer  accurate  because  harmonic  components  of 
the  oscillation  make  important  contributions  to  the  rms 
velocity.  Therefore  A  A  as  developed  in  eq  14  is  typi¬ 
cally  too  small  at  high  q^  values.  This  should  be 
compared  to  the  results  of  Louris  et  al.  [19],  which 
indicate  resolution  decreases  for  q,  values  above  0.73. 

Second,  the  use  of  eq  1  implies  a  more  basic  as¬ 
sumption;  this  is,  the  amplitude  distribution  is  un¬ 
changed  during  the  scanning  of  the  mass  peak.  At  high 
scan  speeds  (i.e.,  >  1)  this  is  clearly  true,  but  at 

low  scan  speeds,  additional  questions  are  raised.  In 
effect,  ions  may  ''diffuse"  from  low  to  high  amplitude 
regions  of  the  distribution  because  of  collisions  with 
neutrals.  This  diffusion  occurs  at  a  rate  that  is  propor¬ 
tional  to  the  reduced  collision  frequency  times  the 
width  of  the  distribution.  Therefore  as  long  as  eq  48  is 
obeyed, 

dA 

—  >  c  A  A  (48) 

at 

the  distribution  will  remain  largely  xmchanged  during 
scanning.  It  is  interesting  to  note  that  by  substitution 
this  expression  can  be  shown  to  be  equivalent  to  the 
condition  >  1. 

At  low-scan  speeds,  this  implies  that  the  ions  may 
actually  "diffuse"  out  of  the  trap  prior  to  the  normal 
ejection  point  of  the  distribution.  In  other  words,  most 
ions  will  enter  a  sufficiently  high  trajectory  for  ejection 
prior  to  the  time  when  the  bulk  of  the  unaltered 
distribution  would  have  reached  the  edge  of  the  trap 


aperture.  This  effect  may  alter  the  ejection  time  for 
ions,  but  it  does  not  alter  the  peak  width  because  the 
relaxation  time  r  is  proportional  to  AA/(5A/50  in 
the  pseudopotential-well  region  and  in  a  diffusive  ejec¬ 
tion  process,  the  peak  width  is  also  proportional  to  this 
relaxation  time. 

At  high  q^  values  the  rate  at  which  ions  may 
change  amplitude  trajectories  is  altered  by  a  secondary 
effect.  The  oscillation  amplitude  now  is  based  on  a 
superposition  of  multiple  frequencies.  Without  colli¬ 
sions  the  phase  relationship  between  the  secular  and 
trapping  frequencies  is  determined  by  the  voltages 
driving  the  oscillations  and  the  initial  conditions,  but 
the  presence  of  collisions  may  alter  the  phase  relation¬ 
ship  between  the  two  oscillation  frequencies.  If  the 
different  oscillations  have  amplitudes  on  the  same 
order  of  magnitude,  it  is  clear  that  minor  changes  in 
phase  due  to  collisions  may  produce  alterations  in 
amplitude  on  the  order  of  A  A  from  a  single  collision. 
In  other  words,  the  collision  frequency  of  the  ions  with 
individual  neutrals  may  determine  the  rate  at  which 
-  ions  diffuse  out  of  the  trap  at  lower  scan  speeds. 
Therefore  the  actual  collision  frequency  c'  described 

by 

m  +  M 

c'  =  nv^a-  =  — — — c  (49) 

determines  the  relaxation  time  constant  r',  and  the 
temporal  line  width  A  t  may  be  expressed  as  the 
smaller  of  eq  20  and  r'  =  4c'. 

These  additional  caveats  are  added  to  the  previous 
geometrical  considerations.  Figure  6  illustrates  a  com¬ 
parison  between  theoretical  predictions  and  recent  ex¬ 
perimental  results  published  by  the  Finnigan  group 
(Louris  et  al.  [19])  and  by  Loundry  et  al.  [3].  Both  sets 
of  results  illustrate  resolutions  between  2  and  4  orders 
of  magnitude  greater  than  the  frequency  line  width 
model  predicts.  Figure  6a  compares  the  maximum 
resolution  prediction  of  the  present  model  to  values 
obtained  by  Louris  et  al.  [19]  for  the  m/z  129  ion  of 
xenon  ejected  at  q^  =  0.73.  It  already  has  been  noted 
that  a  voltage  prediction  of  4.2  V  to  obtain  maximum 
resolution  agrees  very  well  with  the  value  of  5  V 
obtained  by  Louris  et  al.  [19]  at  a  scan  rate  of  5500  u/s. 
It  should  be  noted  that  eq  41  predicts  that  maximum 
resolution  will  be  obtained  at  values  beyond  the 
capability  of  the  commercial  instrument  for  higher 
mass  ions  {m/z  >  100).  Resolution  improvements  for 
high  mass  ions  produced  by  using  axial  excitation 
voltages  between  6  and  25  V  were  recently  demon¬ 
strated  in  our  laboratory  [26], 

Figure  6b  compares  predictions  of  the  present  model 
with  results  recently  reported  by  Loundry  et  al.  [3]. 
The  lack  of  dependence  (which  is  indicated  by  the 
solid  line),  the  lack  of  mass  dependence,  and  the  fact 
that  that  present  theoretical  construction  still  underes¬ 
timates  the  resolution  line  by  a  factor  of  10  are  of 
considerable  concern.  The  authors  believe  this  is  due  to 
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Figure  6.  Comparison  of  the  present  model  with  results  ob¬ 
tained  by  (a)  Louris  et  al,  (19)  and  by  (b)  Lotmdry  et  al.  (3).  The 
solid  line  indicates  predicted  maximum  resolution.  Additional 
details  are  in  the  text. 

a  collection  of  the  previously  mentioned  secondary 
effects.  The  dotted  line  in  this  figure  corresponds  to 

m  (oc' 


This  is  the  expected  scan-rate  dependence  when  the 
temporal  line  width  is  determined  by  the  collision 
relaxation  time  rMn  this  case,  we  took  the  pressure  to 
be  10”“^  torr  and  temperature  to  be  300  K,  The  highest 
resolution  is  obtained  at  =  1  in  this  curve,  and 
the  peak  width  is  determined  by  the  r'  until  eq  20 
again  becomes  the  limiting  expression  for  line  width. 
At  the  higher  scan  speeds,  the  dependence  re¬ 
turns,  but  the  resolution  is  less  than  predicted  because 
of  (at  least)  two  factors  we  have  already  mentioned:  (1) 
the  geometrical  considerations  in  aperture  shape  and 
dipole  excitation  and  (2)  limitation  of  excitation  ampli¬ 
tudes  to  6  V  on  commercial  instruments.  Two  reasons 
for  not  observing  this  same  phenomenon  in  Figure  6a 
also  should  be  noted:  (1)  the  lower  value  in  the 
Louris  et  al.  [19]  data  reduces  the  drive  frequency 
amplitude  with  respect  to  secular  amplitude  and  (2) 


the  initial  amplitude  dispersion  is  wider  for  the  lower 
mass  ions. 

Conclusion 

The  present  amplitude-dispersion-based  additions  to 
the  frequency  line-width-based  model  of  Goeringer  et 
al.  [10]  offer  considerable  insight  into  the  phenomena 
that  affect  resolution  in  ion-trap-based  instruments. 
Furthermore,  it  appears  to  do  so  at  high  scan  speeds 
and  q,  values  that  are  beyond  the  limiting  cases  used 
in  the  construction  of  the  theory,  even  though  consid¬ 
erable  effort  may  be  required  to  develop  the  appropri¬ 
ate  correction  values  required  to  make  the  theory 
quantitative  for  predictions  in  these  regions.  Ongoing 
efforts  by  other  researchers  in  the  areas  of  temperature 
estimation  [21-23],  nonlinear  resonances  [9,  24],  and 
ion  trajectory  studies  [33,  34]  as  well  as  additional 
efforts  in  the  areas  of  geometrical  issues  [27,  28]  and 
collisional  effects  [33,  34]  on  ion  packet  trajectories  will 
all  contribute  to  the  ultimate  understanding  of  ion  trap 
resolution  properties. 
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